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B 

We review recent theoretical progress on glassy dynamics, with special emphasis on the importance 
fj ' and universality of the aging regime, which is relevant to many experimental situations. The three 

main subjects which we address are: (i) Phenomenological models of aging (coarsening, trap mod- 
els), (ii) Analytical results for the low-temperature dynamics of mean-field models (corresponding 
to the mode-coupling equations); and (iii) Simple non-disordered models with glassy dynamics. 
We discuss the interrelation between these approaches, and also with previous work in the field. 
Several open problems are underlined - in particular the precise relation between mean-field like 
^~~- ) (or mode-coupling) descriptions and finite dimensional problems. 
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1 Introduction 

Glassy systems are characterized by the fact that their relaxation time becomes exceedingly 
long for low temperatures, so long that these systems are never in equilibrium on laboratory 
(or even geological) time scales. Notwithstanding, most theories of spin glasses and other 
disordered systems have first aimed at describing the putative equilibrium state of these 
systems, a a In doing so, many difficulties and surprises have emerged - most notably the 
intricate structure of the Parisi solution of the Sherrington-Kirkpatrick (sk) @ mean-field 
model for spin- glasses 13. Despite some early attempts Era, phenomenological and analytical 
descriptions of the non-equilibrium phenomena in disordered systems have only recently 
appeared, which we shall review below. These dynamical approaches have been developed 
mostly because of the accumulating body of experimental data on agingiQ~v3, which is a 
striking experimental consequence of the fact that these systems are out of equilibrium even 
on macroscopic time scales. This aging regime is not the most general out of equilibrium 
situation: a certain degree of universality emerges in the non-equilibrium properties. While 
usual equilibrium dynamics is stationary, i.e. invariant under time translations, the aging 
regime presents a kind of 'covariance': after transients have decayed, the dynamical evo- 
lution of an old system of age t w is described by the same equations as that of a younger 
system of age t w /2, up to a rescaling of time. 

The change of focus from equilibrium to non equilibrium situations also has the in- 
teresting consequence of unveiling strong analogies between disordered systems such as 
spin-glasses, and other types of glasses where disorder is a priori absent, such as fragile 
glasses. These analogies are both phenomenological and formal: many experimental facts 
are in close correspondanc e [ for example, aging phenomena were first studied in detail by 
Struick on polymer glasses 113), but also, the structure of the mean-field equations used to 
describe non-equilibrium spi n-gl asses are almost identical to the Mode Coupling Theory 
(mct) of supercooled liquids l 16 H 17 l. There is thus a strong feeling that the two types of 
systems should be deeply connectedlla ~E9, and there have been several attempts in the past 
few years to establish some precise bridges, which we shall review in this paper. 

The scope of the present paper is mainly descriptive: we focus on general ideas and 
concepts rather than on more technical aspects. We refer the reader to the relevant papers 
for more details. 

1.1 Experiments: History dependence and Aging 

The simplest way to see that spin-glasses below the phase transition temperature T g are 
not in equilibrium even after times of the order of hours (or more) is the following: the 
sample is quenched rapidly (under zero magnetic field) from high temperatures T 3> T g to 
the working temperature T\ < T g which is reached, by convention, at time t = 0. Then a 
very small oscillating field is applied to measure the a.c. susceptibility \ °f the sample at 
a certain frequency uj. What is observed is a slow continuous decrease of the amplitude of 
X as a function of the time L^ elapsed since the sample reached the temperature T±, which 
is called an aging effect. OtISI In other words, x is a function of both frequency and time: 
x(w, t w ). The response of the system to a perturbation thus depends on the thermal history. 
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Figure 1: The amplitude of the out of phase magnetic susceptibility x"(w,t w ) vs cut w for the insulating 

spin-glass CdCr1.7Ino.3S4 (from Ref. [14]). The frequencies are u) — 0.01, 0.03, 0.1, 1. Hz and t w is the total 

time elapsed after the quench. In the inset, a more refined scaling for x"(^,t w ) as discussed in Ref. [14]. 



To a good approximation, the shape of x(w,i w ) can be parametrized as follows (see 
Fig.l): 

X(u,t w ) = A (ujt w )- b + X st(u) (wf M > 1) (1) 

where A is a temperature dependent amplitude and b an exponent that moves in the range 
0.1 — > 0.4. The important points of the above parametrization are: 

• The response function is the sum of a stationary part Xst(^) which is independent 
of the age of the system t w , and of an aging (or non-stationary) contribution, which 
decreases with time. Xst(^) behaves as oo a with a small exponent a (sometimes called 



, or perhaps as logo;. For systems in equilibrium, the time dependent (aging) 
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contribution disappears. 

• The aging contribution can be described with a function of the scaling variable l 13 H 14 l 
Lot w . In general, the susceptibility of a system with a single relaxation time r is a 
function of lot. Hence, the above scaling form means that the effective relaxation time 
of the system is of the order of its age t w itself. (SeeE-3 for a discussion of the inset in 
Fig. 1 and a more detailed discussion of an alternative description of the Xag data.) 

Another set of experiments which basically carry the same information is those of the 
so-called 'Thermo-Remanent Magnetisation' (TRM) relaxation l 13 H 14 l. The system is cooled 
under a small magnetic field H , which is left from t = (the time of the quench) to t = t w , 
and then suddenlvswitched off. The subsequent relaxation of the magnetisation M can be 
decomposed asEjIlj 



M(t w + r, t w ) = M st (t) + M AG (t w + r, U 



M 



st (t) = lim M(t w +r,t v 

t-w—tOO 



(2) 
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Figure 2: The aging part of the thermoremanent magnetization M& a (t w + T,t w ) (normalized by the zero 
field cooled value M{ c ) vs. r/t w for AgMri2.e, in a log 10 scale (from Ref. [14]). The sample was cooled in a 
0.1 Oe field from above the glass transition T a — 10.4K to a subcritical temperature T — 9K. It waited for 
t w — 300, 1000, 3000, 10000, 30000 sec. under the field that was suddenly switched off at t w . The decaying 
magnetization was recorded during all subsequent times t + t w . 
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where again, there is a 'fast' ]] stationary contribution M st (t) which is independent of t w , 
and an aging part M AG (t w + r, t w ) which to a good approximation (see Fig. g and Section H() 
is a function of the ratio t/(t + t w ). This again suggests that the effective relaxation time 
of the system is of the order of its age. Actually, within the linear response approximation, 
x(uj,t w ) and M(t w + r,t w ) are essentially Fourier transform of each other. More precisely, 
introducing the response function R(t,t'), one has 



M(t w +r,t w ) = H dt' R{t w +T,t' 



x{UiU 



dt' R(t w ,t')e 



itu(t'—t w ) 



(3) 
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TRM relaxation in spin-glasses t3 and stress relaxation I13B3, electric polarisation ElEJ or 
specific heato relaxation in many very different glassy materials show - rather remarkably 

- similar features, with a fast initial drop at small times r, followed by a slow decrease of 
the signal on time scales of the order of the waiting time t w . The same picture also pertains 

- on much smaller time scales - to numerical simulations of the response function of the 
three- c3n3 and four- dimensional c3 Edwards- Anderson and in mean-field C3~V3 spin-glass 
models. 

Since the response function depends on some aspects of the thermal history, it is inter- 
esting to consider more complicated experimental pjxiiocoles such as different cooling rates 
from high temperatures to Ti, temperature cycles E3H between two temperatures T\ and 
T2, or even field-cycling Ei. The detailed discussion of these situations is beyond the scope 
of the present paper. However, it is interesting to notice that while the a.c. susceptibility 



It is fast in the sense that experimentally, the major part of M st (t) has decayed to zero after the first 
second. However, M S t(t) only decays as (to/r) a with a small: see after Eq. (hi) above. 



x(w, t w ) depends extremely weakly on the cooling rate in spin-glasses, there are experimen- 
tal systems (e.g. dipolar glasses EJ) for which this dependence is large. We shall come back 



to this point in Section 2.2 



Finally, let us mention that aging can also be seen in correlation functions (rather than 
response functions). In equilibrium, these two quantities are related by the well-known 



Fluctuation-Dissipation Theorem (fdt), which is, as we shall discuss in Section 1.2, not 
necessarily valid in out of equilibrium situations: in general, correlation and response do 
not contain the same information. 

From an experimental point of view, correlations are obtained from time dependent 
noise spectrum measurements E3'E3 S(io,t w ), which are rather more difficult than response 
measurements fl From a numerical point of view, however, it is very easy to compute 
the time dependent correlation function C{t w + T,t w ), which for a spin system is defined 
as C(t w + r, t w ) = jjJ2iSi(t w + T)Si(t w ). The behaviour of C(t w _+- r,t w ), obtained on 
relatively short time scales from simulations for finite dimensional c3c3 and mean-field£j"EH 
spin-glass models, reveals aging below the spin-glass transition (for reviews, seeE3c3), again 
qualitatively described by a shape similar to Eq. (||), with M's replaced by C's. 

1.2 Think in the two times plane : time sectors. 

Out of equilibrium situations have long been considered as untrustworthy. What the exper- 
iments tell us (and what the theories below will confirm) is that provided one abandons the 
idea that the correlation or response functions should be 'time translational invariant' (tti), 
one can make sense of the experimental data by explicitly keeping the dependence on the 
two times: C(t w + r, t w ) ^ C(t). Equivalently, the Fourier transform of these quantitites 
will not be functions of the frequency oo only, but of both uj and the time since the quench 
t w . 

More precisely, for any physical system there are a priori two other time scales, one of 
them is microscopic (and will be noted t ) and determines for example, the single spin-flip 
time. The other time scale is the equilibration time £ E rg> which, for a finite size system, will 
always be finite (albeit often astronomical). The regime in which one expects to observe 
some 'universal' features (independent, for example, of the details of microscopic dynamics) 
is the regime where: 

to < t w + r < t ERG and t < t w < * E rg • (4) 

This does not require any particular relation between r and t w . For t ERG = oo, one may in 
general expect (and one indeed finds in some models) that the functional form of - say - 
the correlation depends on the way t w + t and t w are taken to infinity. 

The simplest example to see this is the case of ferromagnetic domain growth. After a 
quench from high temperatures to a non-zero temperature below the Curie point, a pattern 
of domains of positive and negative magnetizations starts coarsening. The typical domain 
size £(t w ) at time t w diverges as a power law (or possibly more slowly in the presence of 



c This requires to perform many independent quenches where the magnetic noise is recorded for different 
ages tw and then averaged over the different quenches. 
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impurities - see Section p^). The ergodic time t ERG is then the time at which the size of 
the domains is that of the sample. We shall consider the thermodynamic limit in which 
£erg = oo. Within an approximate (large n) theory of coarsening, the correlation function 



for large times is indeed found to beta 

C(t w + r, t w ) = C st (t) + C AG ( J} tw ) r ) ■ (5) 

The first term describes the fast relaxation of the spins within each domain, and has 
the same form as it would have in equilibrium, when there is only one infinite domain. Its 
limitf]: 

q EA = lim lim C(t w +T,t w ) (6) 

is an example for this simple model of the Edwards- Anderson parameter, a quantity that 
plays an important role in glassy dynamics (in this case it is simply the magnetization 
squared). The second term of Eq. (J5|) describes the relaxation of the system due to the 
motion of domain walls and it manifestly depends on the waiting time. 

As t w goes to infinity, one will be probing two distinct regimes, depending on whether one 
takes t w — > oo with r finite ('stationary' regime) or T,t w — > oo with 
£,(tw)/£(tw + t) < 1 ('aging' or coarsening regime). In terms of the correlation function, the 
stationary and aging regimes are simply defined as the regimes of (large) times in which 
C(t w + r, t w ) > q EA and C(t w + r, t w ) < q EA , respectively. 

These considerations can be translated in Fourier space as follows. Defining 

C(w, t w ) = / dr C(t w ,t w - t) e WT , (7) 

Jo 

one finds that C(cu, t w ) does not only depend upon to (as it would in equilibrium situations) 
but also upon t w . For example, for the correlation of the form (||), one obtains, in the limit 
urt w » 1: 

<>(«,*.) -<VM + |Kw) <8) 

where C is a certain function related to C AG . In the simplest case in which £ grows as a 
power law, one finds £(t w )u / £' (t w ) = u>t w . 

More generally, one can envisage the possibility that different physical mechanisms act 
on different large-time sectors, defined as 

f = 0(l), r^T = 0(D, ^=0(1), etc. (9) 

to hi(t w + T) h 2 (t w + T) 

where the different functions hi (no longer necessarily related to domain sizes) are monotonously 
increasing functions which grow differently, in such a way that: 

< MU < 1 ^ hj(t w ) =1 for . < 

hi{t w + T) hj(t w + T) 



^Note the crucial ordering of the two limits. 



Notice that these time sectors correspond to asymptotically distinct relative 'epochs' in the 
sense that if t\ , ti belong to the domain defined by hi and £2 , *3 to the one defined by hj 
with j > i, then ti, £3 also belong to the sector by hj. 

The correlation function can be, for example, the sum of terms of the type (||), each 
with a different scaling function hi replacing £. Because the different scalings vary in time- 
sectors that do not overlap, such a function cannot be reexpressed in terms of a simpler 
scaling form valid for all large times. 

A simple example for the hi(t) is 



(1 - m)to~ 



Ht) = exp ( -x^- (11) 



with < m < 1. In this case the i scaling form corresponds to r ~ t^^tlfi. In particular, 
H = yields the time-translational invariant form and fi = 1 the r/t w scaling variable which 
is independent of the microscopic time to- 

The main message of this Section is that once one abandons time-translational invari- 
ance, as one should in systems that never equilibrate, the two-time correlation (or response) 
function in the long-times limit may have a rich structure including multiple-scaling forms 
like (J9|) . It will turn out that a rather general classification of the asymptotic behaviour in 
the two-time plane can be made using only the monotonicity property of the correlation 



and simple group theory (see Section 3.4). 

The simple example of coarsening also illustrates the fact that in order to decide whether 
a system can be considered to be in equilibrium, one-time quantities (such as the energy, 
magnetization, etc) can be misleading. Indeed, at any finite time the excess energy density of 
a coarsening ferromagnet is proportional to the total domain surface divided by the volume, 
a quantity which soon becomes very small. If one were to judge the degree of equilibration 
by only measuring the excess energy density, one would wrongly conclude that the system 
equilibrates rather rapidly. On the other hand, two time quantities (such as the correlation 
function) reveal very clearly that the system is still out of equilibrium even at long times. 

1.3 Fluctuation- Dissipation Relations 

As was mentioned before, in a system at equilibrium, the response to an external magnetic 
field and the autocorrelation functions are related through the fluctuation-dissipation theo- 
rem (fdt). This is true in general for the response to a field h conjugate to any observable 
O and the corresponding autocorrelation C (t w +T,t w ) = (0(t w -\-r)0(t w )} . In equilibrium: 



tto{tw + T , tyj) = 



1 <9C , E q(t 

= ^0,EQl"TJ = -■ 
h=0 



iWr) = -- ^ w (12) 



5h(t w ) 
or, introducing the integrated response Xo{t w + t, t w ) = J t w T Ro(t w + T, t')dt': 

Xoit w + r, t w ) = XO,eq(t) = — (Co,eq(0) - Co,eq(t)) . (13) 



If O is the energy, one obtains the relation between energy fluctuations and specific heat, 
if O is the magnetization, one finds a relation between the (time dependent) field induced 
magnetisation and the magnetic noise correlations, etc. 

In order to study the relation between 'fluctuations' and 'dissipation' in out-of-equilibrium 
systems, one has to think in terms of two-time correlation and response functions. Let us 
then consider a given t w and make a parametric plot of Xo( T + t w ,t w ) vs. C (r + t w ,t w ) 
when r varies. One then takes a larger t w and repeats the plot, and so on. If the system 
equilibrates after a finite time £ E rg! one obtains, when t w 3> i BRG , a limiting Xo vs. C curve 
which is a straight line with slope — 1/T: this is the fdt. 

Consider instead what happens in the example of domain growth in an infinite size 
system. Within a large n treatment of the problem e3, one obtains for the x vs - C curves a 
family of curves shown in Fig. |3|. For large t w , the curves approach a broken line: one with 
slope — 1/T for values of the correlation larger than q EA (i.e. for times in the stationary 
regime defined in the previous Section), and one with zero slope for values of the correlation 
smaller than q EA (the aging regime). This simple example illustrates that (despite the fact 
that the dynamics is very slow when t w — > oo) the system cannot be thought of as in a 
'quasi-equilibrium' state, for which concepts from equilibrium are more or less valid: there 
is always a regime in which fdt is strongly violated. (It would be interesting to confirm 
these results within more realistic models of coarsening, and also, obviously, experimentally 
and numerically.) We also see that the Edwards- Anderson parameter plays a role for the 
response functions. 

Similar results are found analytically in the mean-field spin-glass systems which we shall 
review below (see Fig. H and Section |3|). The main difference is that the aging part of the 
curve (i.e. C < q EA ) has a non-zero slopeci~E3 This is quite important, since it means that 
the integral of the response function over the aging regime gives a non-zero contribution. 

In the more realistic 3D Edwards-Anderson model, the form of the x vs - C curves 
which are obtained numerically (at least for the computer times accessible at present) are 
actuallv remarkably similarE3 to the mean- field prediction for the corresponding mean- field 
model e3. The same is true of some very recent numerical simulations of binary soft-sphere 
mixtures c3. 

The violation of fdt can be parametrized by introducing a violation factor X (t,t') 
defined as 

Mti0a *.ftOW) , (14) 

Note that we are differentiating with respect to the smallest time t' . In analytic studies of 
mean-field systems, one can furthermore show that for large times X depends on t, t' only 
through the value of the correlation function: X (t,t') = X[C (t,t')}. In particular, when 
Co > Qea,0) ^o = 1, and the FDT is recovered. 
It turns outc3 that the 'effective temperature' 

T EFF (^ + r, t w ) = T (15) 

^-cnf-ui ~r T, l w ) 

is precisely the temperature which would be read on a thermometer with response time r (or 
frequency u ~ 1/t) when connected to the observable O at time t w . A 'fast' thermometer 
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Figure 3: The susceptibility x(£u> + t, t w ) = J w T ds R(t w +r, s) vs the correlation C(t w + r, t w ) for domain 
growth in a n-componcnt vector ferromagnet at T < T a , obtained from an analytical treatment of the large 

n limit. (From Ref.[48].) 

of response time r <C t w will then probe the stationary regime for which X = 1 and thus 
measure the heat-bath temperature. This is the reason why glasses, although still out of 
equilibrium after many hours, feel as cold as the room they are in. 

1.4 Edwards-Anderson parameter, weak-ergodicity breaking and donation 

The situation we have described in the last two subsections is one in which equilibrium is 
not achieved, in the sense that configurations are not visited with a probability given by the 
Gibbs-Boltzmann weight. However, this is not the main point: none of the above results 
can be explained within a strong ergodicity breaking scenario where the system falls into a 
very long-lived metastable state, and achieves fast equilibration within such restricted sector 
of phase-space. A mixture of hydrogen and oxygen, or a crystal of diamond are systems 
which are metastable, but for all dynamical purposes in equilibrium. As Feynman puts it: 
'Equilibrium is when all fast things have happened, and slow things not yet'. The out of 
equilibrium situations which are of interest to us are those where, in a sense, 'things keep 
happening on all time scales'. 

More precisely, in a spin-glass (or a coarsening problem) in the low temperature phase 
the spins do, on average, remember for some time their orientation, which leads to a non- 
zero Edwards- Anderson parameter q EA = lim T ^ 00 hm tiu ^ 00 C(t w + T,t w ). However, if the 
waiting time is finite, the system is able to escape arbitrarily far from the configuration it 
had reached at t w , leading to 



lim lim C(t w + r, t w ) = 

t w — >OOT—KX) 
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(16) 
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Figure 4: The susceptibility % vs the correlation C for the 'discontinuous' mean-field spin glass model of 

Section g (p = 3) at T < T B . In this case the fdt violation for C < q EA ~ 0.76 is given by —X/T, with 

X = (I — qEh)/<iEA ~ 0.76. More complicated situations are described in Section pi (From Ref.[48].) 



even in the low temperature phase where q EA > 0. This situation was called 'weak-ergodicity 
breaking' in Refs.p3'El'E3]. 

It is important to note that the fact that a system undergoes weak (as opposed to 
strong) ergodicity breaking for infinite times does not mean that stable states do not exist 
in phase-space. The dynamical behaviour depends on the choice of initial conditions. The 
example of ferromagnetic coarsening is again eloquent in this respect. There exist two true 
equilibrium states, and the relaxation from an initial condition close to one of them (e.g. all 
spins up) is fast and shows no aging. But in the relaxation from a random initial condition 
(e.g. a typical configuration at a temperature above T c ), the system remains forever (in the 
thermodynamic limit) undecided as to which state it will go to. 

Another interesting question is the following: suppose that at time t w one 'duplicates' 
the spin system and evolves subsequently the two copies using two independent thermal 
baths. (This process was called 'donation' in Ref.pl].) Will the two copies stay 'close 
together' or conversely will they evolve independently, forgetting their common breed? This 
is measured by the 'overlap' function Q(t w + r,t w + r) defined asH 



N 



Q(t w + r,t w + r) = -Y^ Sl L) (t w + r)S\ l > (t w + r) 



:(2), 



iV 



sl l \t w ) 



sf\t w ) 



(17) 



where the superscripts ( 1,2 > refer to the two copies. We shall see below that the limit 



lim lim Q(t w +T,t w +r) = Q c 

t w — >00 T— »00 



(18) 
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can be zero or non zero ev en if » E » is non zero. This may serve to distinguish different 
types of aging dynamics EJ'eSH 30 !! 31 !. 

If the limit is taken in reversed order (i.e., t w — ► oo first), thejiverlap function contains 
the same information as the correlation function. More preciselyt^: 



lim Q(t w + T,t w + r) = C st (2t) . (19) 



2 Phenomenological models of aging 

In order to account for the above experimental results, one has as usual the choice between 
some physically motivated, but phenomenological, pictures or some rather more precise mi- 
croscopic models, in a limit in which they are analytically tractable. Both approaches are 
actually complementary, and shed light on each others' limitations. We shall start by re- 
viewing the phenomenological pictures of aging, based either on domain growth arguments, 
or on models of random walks in phase space. 

2.1 Coarsening in non disordered systems 

As we noted in the previous Section, the simplest model where aging occurs is the ferro- 
magnetic Ising model suddenly quenched below its Curie point temperature. The initial 
configuration is random, and it orders progressively through domain growth. Depend- 
ing on the class of microscopic dynamics t3, the typical size £ of the domains grows as 
t ' 2 ('non-conserved' case) or t ' 3 ('conserved' case). The 'age' of the system is thus 
directly encoded in the spatial correlation functions. The two-time correlation function 
C{t w + T,t w ) = jjJ2iSi(tw + T~)Si(t w ) can be calculated exactly in some cases (e.g. the 
Ising model in one dimension, or the large-n 'spherical' model), or using some approxima- 
tions (for a review see Ref-pHI]). One finds an expression as Eq.(||) 

C(t w + r, t w ) = C st (t) + C AG ( J} tw ) r ) (20) 

where C AG decays as a (non trivial) power-law for large arguments. C AG (1 -Liz) is in general 
also singular around u = 0, its behaviour is characterised by an exponent e3 b and, in the 
'non-conserved case', one can argue thatE3S, C AG (1) — C AG (1 + u) oc y/u. Correspondingly, 
the aging part of the a.c. susceptibility decays as in Eq. (|l]), with b = 1/2. 

The overlap function Q(t w +T,t w + r) can also be estimated in this simple coarsening 
situation. Within standard approximate treatments of coarsening eSpEH, one finds that the 
quantity Qoo defined in Eq.fllq) is non zero, meaning that the two copies follow each others' 
footsteps (within finite times) in their evolution towards equilibrium. Note also that within 
the large n approximation, the fdt violating factor X goes to zero at infinite times, as 
£ _1 . It would be interesting to know whether this is a more general property of coarsening 
dynamics. 

It is interesting to remark that in the presence of a small external magnetic field H, 
aging is 'interrupted' after a finite time t ERG (H), since one of the two phases is favoured 
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by the field. For example, in the case of the spherical model, one findsa t ERG (H) oc H~ 2 . 
This behaviour is expected on general grounds: it corresponds to the time beyond which 
the curvature induced driving force is superseded by the field induced driving force. 

2.2 Coarsening in disordered systems I: Random ferromagnet or random field 

Let us now consider the case of a disordered ferromagnet in dimension larger than 2, where 
random local magnetic fields or random local couplings are present. If the disorder is suf- 
ficiently weak (for example, if a small fraction of the ferromagnetic bonds are removed), 
the ground state of the system still has long-range ferromagnetic order, and the description 
in terms of the growth of ordered domains is valid. However, due to the presence of dis- 
order, the domain walls will tend to be pinned by local inhomogeneities. The problem of 
domain walls in disordered environments has been the focus of intense study in the recent 
years E3. In many aspects, this problem is close to the spin-glass problem, with a large 
number of metastable states (although of course the 'spin-glass' nature of the problem only 
concerns the small fraction of spins which belong to the domain walls). The dynamics of 
a given section of a domain wall proceeds by thermally activated hops between different 
favourable configurations. Ordered domains thus grow on average, but at a much reduced 
rate compared to the pure ferromagnet described above. 

A generally accepted description is as follows E3: the typical pinning energy scale of 
a domain wall of linear size R grows as TR e , where 9 is an exponent which depends on 
the problem (random field/random bond) and the dimension of space. A simple scaling 
argument! then suggests that the time needed for a domain to reach a certain size R is 
given by 

r(H)oct exp(— J . (21) 

After a certain time £, the typical size of the domains is thus expected to be given by 

i 
' kT log (fY 
«*)«| 1^| , (22) 

provided the corresponding pinning energy T£(£) is large compared to kT. (In the other 
limit, the pinning energy is negligible, and one recovers the growth law which we discussed in 
the previous paragranh.) This logarithmic growth of domain sizes has been rather carefully 
checked numerically Ej in D = 3 (T / 0). 

Apart from the fact that domains grow very slowly, one expects that the picture prevail- 
ing in the pure ferromagnetic case is not drastically modified. In other words, the correlation 
function should still be given by Eq. (0), but with a logarithmically growing £(£). This 
was confirmed numerically on the Ising model with random fields, for D = lS and D = 3 
l3. In D = 1, the motion of a 'domain wall' (a point) is given by Sinai's diffusion law, i.e. 



e This argument assumes that barrier hetsbte between metastable states behave (as a function of R) 
similarly to the pinning energy of each state t_ll£J 
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£(£) oc log (t) (at least for times smaller than a certain temperature dependent t ERG , which 
diverges when T — > 0). The aging part of the two-time correlation function can indeed be 
satisfactorily rescaled when plotted versus £(t w )/£(t w + r). 

It is important to notice that the above scenario, where the system tries to reach a well 
defined (temperature independent) state, but is slowed down due to pinning by impurities, 
leads to large cooling rate effects. This is because the crossover energy T£(£) ~ kT will 
be reached at later times if the cooling rate is slower Ej. The size of the domains when the 
working temperature T± is reached will thus be larger, the smaller the cooling rate, and this 
will affect many physical observables, such as the energy. 

Let us finally point out that disorder is actually not necessary to obtain logarithmic in 
time (rather than power-law) growth of the domain size. This was first shown in Ref.pl] for 
a pure Ising model with next-nearest neighbour couplings. Actually, if the domain walls are 
below their roughening temperature, the dynamics proceeds via the nucleation of terrace£3, 
and also leads to logarithmic domain growth. 

2.3 Coarsening in disordered systems II: Spin glasses and droplets 

It is not obvious whether the above coarsening description also applies to spin glasses, 
because of the non-conventional nature of their order parameter. However, Fisher and Huse 
Ej have argued that for spin-glasses in finite dimension and for any given temperature below 
the spin-glass transition, there are only two 'pure states' (spin-reversed from each other) 
which have to be considered, which can conventionally be called 'up' and 'down', very 
much like in the Mattis modela. This assumption g ives s ome physical content to a scaling 
description of the spin glass phase, first advocated inrir 6 !. It allowed several authorsmrH 68 ! 
to develop a rather complete phenomenological picture of spin-glasses in low dimensions, 
where the spin-glass is considered as a 'disguised ferromagnet' (with however the important 
difference that the two pure states are not stable when the temperature is changed -,-see 
below). This is at variance with the mean-field picture emerging from Parisi's solution of 
the Sherrington-Kirkpatrick modelEI, where many (non trivially related) pure states coexist 



If this 'two-state' picture is retained, the dynamics of the system can again be described 
in terms of growing and coalescing compact domains^, which have also been called 'droplets' 
l3 in this context. The presence of disorder presumably pins the domain walls, leading again 
to a logarithmic growth of the droplets, and thus to a two-time correlation very similar to 
the random-bond or random field case described above. In particular, the aging part of the 
correlation function should be a function of \\ £(t w )/£(t w + r). 

One striking property of experimental spin-glasses, however, is the very weak depen- 
dence of its physical properties on the cooling rateE9. For example, the asymptotic value 
of the a.c. susceptibility, x(o;, t w — » oo), is nearly independent of the cooling rate. This is 
at first sight surprising in a scenario of activated domain gr owth, as we emphasized above. 
But if one argues that the spin-glass phase is 'chaotic' r 6 H 63 H 67 l, i.e. that the two pure states 
towards which the system evolves are extremely fragile to temperature changes, it becomes 



•^ Fisher and Huse actually postulated that the scaling variable would rather read £(t)/£(£™). 
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obvious that the dynamics of the system at temperatures greater than T\ is useless to bring 
the system closer to its equilibrium at the working temperature T\. In a first approxima- 
tion, the configuration reached by the system at temperature T\ + 8T is as remote from the 
'true' equilibrium state at temperature T\ than a high temperature configuration. Hence 
the cooling rate may indeed have a negligible effect. 

How does the two-state picture compare with aging experiments or numerical simula- 
tions, in the case of a 3 — d Ising spin-glass? The detailed discussion of this po int is beyond 
the scope of the present paper and the conclusionsstill controversial aOOtla, in particular 
the existence of a transition in a magnetic field L3~Ea. One should however mention the 
following results: 

Both experimentally and numerically, the aging part of the correlation and re- 
sponse functions follow a scaling which is systematically closer to r/t w than the expected 
£(t w )/£(t w + r) (or even £(T)/£(t w )) if the barrier heights scaled as R e , e.g. £(t w ) = 
log(t w /t ) I . The point is that such a scaling would place the curves of Fig. 2 in the 
reversed order: the young ones would be below the old ones in a r jt w plot. 

However, as suggested by RiegerE3, the two state picture with barrier heights scaling as 
T log R (i.e. Q_— > 0) would lead to an algeb raic d omain growth law (as postulated by Koper 
and Hilhorst c3) : £(£) oc t a with a = A;T/TE3Ea, and thus in turn to a r/t w scaling. 

A direct numerical indication of a growing length scale was searched for inEaE3. A 

possibility is to study the following correlation function 



G <?>*) = m E^w^w-O^-O*)) > ( 23 ) 



i 



where S\ and S\ which has the following intuitive meaning: knowing that the copies 1 
and 2 are in the same state at site i (resp. opposite states), what is the probability that they 
are still in the same state (resp. still in opposite states) a distance r apart ? Numerically, 
G(f,t) is seen (in 3 dimensions) to be of the formO: 

° m = Vt g {it)) m *e (24) 

which indeed suggests the presence of a growing scale £(£) in the dynamics. The conventional 
droplet picture predicts (apart from a logarithmic, rather than power-law, growth of £(£)) 
that C, = 0, since the equilibrium state should be unique up to a global sign change, whereas 
simulations by Marinari et al. IZ3E3 suggest that £ > 0, just as in the equilibrium situation 
obtained with replica field theory L3 in dimensions smaller than 6. 

The temperature cycling experimentsEj where the system is cooled to temperature 

T\ < T g , then to temperature Ti < T\, and finally back to temperature li, show a very 
striking conjunction of 'rejuvenation' (when the temperature is decreased) and memory 
(when the system is heated back) . Th e coe xistence of these two effects is rather awkward 
to interpret within the droplet picture l 13 H 79 l. In the same spirit, Weissmann et al. E3 have 
argued that the 'second noise spectrum' of spin-glasses does not conform with what could 
be expected from a simple 'two-state' picture. 
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In summary: even if the 'disguised ferromagnet' picture did provide a correct description 
of the equilibrium properties of low-dimensional spin-glasses, it is not obvious that this 
description is sufficient to account for the out- of- equilibrium properties. One reason is 
that - precisely because of the chaotic nature of the equilibrium phases - the system will 
not only nucleate domains of the nominal equilibrium state, but probably also 'phases' 
corresponding to nearby temperatures, which will thus contribute to the non-stationary 
part of the dynamics. 

2.4 Stranded in phase space: the 'trap' model 

'Phase-space' models are another very useful class of phenomenological models for the 
dynamics of complex systems, and have been advocated by very many authors over the 
yearir°lr 5 l. The dynamics of the whole system is summarized in the motion of a single point 
evolving within a complicated energy landscape in configuration space. From a general point 
of view, one expects this energy landscape to be made of 'valleys', or 'traps' (within which 
all configurations are mutually accessible in a short time) separated by 'barriers', which the 
system can only overcome by thermal activation. A coarse-grained representation of the 
problem can thus be given in terms of states a,/3, 7... between which the system wanders. 
The dynamics of the system is thus described by a master equation for the probability to 
find the system in the state a: 

dP 

"of = " £ W <*^ P °> + E Wp.+aPf) • (25) 

The choice of the hopping rates W a -+8 then encodes the statistics of the barrier heights and 
the geometry of the phase space. It is rather arbitrary apart from the constraint of detailed 
balance. For example, one can organize the 'traps' on a hierarchical tree, and choose W a ^p 
to only depend on the distance between a and (5 along the tree. This has led to several 
'ultrametric' diffusion models E3, with many interesting results, including, in some cases, 
aging effects e3~e3 

• The one-level tree. 

The simplest of these modelspa for which the appearance of aging has a particularly 
clear interpretation, is when the hopping rate only depends on the starting state: W a ^p = 
(iW a ) _1 , where N is the total number of states. (The final state j3 is thus independent of the 
initial state; the process starts anew at each jump.) This corresponds to the picture drawn 
on the left of Fig. ||. The trapping times r a are of the form t exp(B a /(kT)), where B a is 
the energy barrier 'surrounding' state a. Within this description, the equilibrium measure 
P^ q (if it exists) is simply proportional to r a . In order to reproduce the correct Bolzmann 
equilibrium, one should thus identify B a with the free-energy of the statef] f a . Mean-field 



9 The relation between barrier heights and energy depths is however not obvious in general; for recent 
work on this subject, see Ref. p] 
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models of spin-glasses EH or replica treatment of randomly pinned mani 
distribution of the metastable states' free energies f a is exponential E9' 



suggest that the 



p(/ a )ocex P (-^) (26) 

with a certain parameter x < 1 in the glassy phase, x = T/T g in the Random Energy 
ModeE3'0. The appearance of this exponential tail for 'deep' states can be understood on 
general grounds, and is related to the so-called 'extreme value statistics' (for a more precise 
discussion, see Refs. pac§|). The corresponding distribution of trapping times |j is then 
easily found to be 

p{r)dT = p(f)df — p(r) oc r » to -^ . (27) 

Let us now introduce the quantity H(t w +T,t w ) defined as the probability that the system 
has not changed trap between time t w and time t w + r. This quantity is found to be 
very different depending on whether x is larger or smaller than 1. In the former case, 
lim^^oo H(t w +T,t w ) is well defined, and found to be proportional to (t /r) x . For x < 1, 
however, one finds that Ii(t w + r,t w ) 'ages', and is given by: 

n(i w + r, t w ) = ^^ f 1 du{l - uf- l u- x (t w » Q . (28) 

IT J-£— 

T+t m 

In physical terms, this means that after a waiting time t w , the only states which have 
an appreciable probability are those with a trapping time of the order of t w itself. This 
reflects the fact that the distribution of trapping times p(r) becomes so broad when x < 
1 (the average trapping time becomes infinite), that the sum of all the trapping events 
T i + T 2 + • • • + T N is always dominated by its largest term, which is thus of the order of the 
experimental time itself E^ .This dominance of a few very important events is a characteristic 
feature of Levy statistics^. On the other hand, when x > 1, the trapping times are all of 
order t . 

Let us now define the overlap q a/ 3 between states as q aa = q EA and q a ^p = qo (more 
general choices will be discussed below). The self-overlap q EA is smaller than 1 in general; 
this reflects the fact that many microscopic configurations are mutually accessible within 
times of order t ; this contributes to the equilibrium part of the correlation function which 
the one-tree level cannot describe. The spin-spin correlation function, averaged over the 
disorder j] is thus given by 

C AG (t w + r, t w ) = q EA U(t w + r, t w ) + q (l - U(t w + t, t w )) . (29) 

Note in particular that lim^oo lim^^oo C AG (t w + r,t w ) = q EA when x < 1, but that 
lim^oo ]imt w -*oo C AG (t w +r, t w ) = qo for x > 1. Within this model, x = 1 thus corresponds 
to a true glass transition. 

h Generalisation to other distributions of f a has been considered in Ref.^ZI] 

1 i.e. over the distribution of f a . Note that in this model, the correlation function is not self-iweraging, 
precisely for the same reason as for the static overlap distribution P(q) in mean-field spin-glasses Ej 
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Figure 5: Schematic phase-space landscape of a one level tree and of a multi-level tree. 

Equation (|2q ) leads to the following asymptotic behaviour for C AG (t w + r,t w )^B: 

sin(7rx)(g EA - q 



^ Q0 + 



7r(l — x) 

sm(Trx)(q EA -q ) ft 



l-x 



TT.r 



r < t w ) (30) 

(r » £„,) (31) 



Hence, both the 'short time' (t <C t <C t w ) and 'long time' regimes are described by 



power-laws, much like in the simple coarsening models described in Sections 2.2 and 2.3 



The above model can be endowed with magnetic properties by assigning to each state 
a certain magnetisation m a , and modi fyin g the hopping rates in the presence of a field to 
recover the correct equilibrium weights E3. One then finds that the generalised form of the 
FDT (Eq. (|14])) holds, with: 



X(t w + T,t w ) = 1-C + C 



t w + T 



(32) 



where C, is a free parameter of the model, restricted to the interval [0,1]. For C, = 0, the 
thermoremanent magnetisation is simply proportional to C(t w + T,t w ). 

Experiments on the thermoremanent magnetisation or the a.c. susceptibility show that 
both the 'short time' and 'long time' regimes can be fitted by power-laws (see Fig. 7 in 
Ref.^5|). The exponent x which comes out of these two fits is however different: x = 1 — b ~ 

0.6 — 0.9 (depending on temperature) from the r<Ct ffi region, and x ~ 0.1 — 0.3 from the 

^ j. • I79II14I 

t ^> t w region liaO. 

One can also define the overlap function Q(t w + r,t w + r) within such a model, and 

show that it is simply related to C through Q(t w + T,t w + r) = C{t w + 2r,t w ). One thus 

finds, in particular, Qoo = qo: two copies of the same system decorrelate completely from 

each othertiil. Note however that this would not be true if the 'traps' were organized along 

a low dimensional 'path' in phase space. 
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• The multi-level tree. 

An interesting generalisation of this 'one-level' trap model is to consider a hierarchical 
organisation of traps within traps, and to relate the overlap q a p between two states to their 
distance along the tree (see the figure on the right of Fig. |5|). The hopping rate W a -tp is 
still taken to be independent of the final state for a given distance along the tree, or for a 
given overlap q = q a $ 

W a ^ oc — . (33) 

As reviewed in Refs pEiJa] , the interpretation of the static full replica symmetry breaking 
solution suggests that the distribution of free-energies of the states at a certain level of the 
tree is still exponential, but with a parameter x which now depends on the overlap between 
these states (and is the inverse of the Parisi function q(x)) 

ftCfl <* exp (-^PQ) , (34) 

where / is the free-energy of the 'ancestor' state, itself distributed exponentially with a 
parameter x(q — dq), etc. Assuming again that the barrier between the states is related to 
the free-energy depth as pictured in Fig. [BL one is led to surmise that the trapping time 
T a , q is still distributed as a power-law, but with a q dependent exponent: 

Ml) 
P*( T ) ^»t ^^ ■ (35) 

Note that x(q) is an increasing function of q, which means that the smaller the overlap 
between states, the broader the distribution of time scales. In other words, 'fast' processes 
are deep down the tree. The correlation function is now determined from| 

M M 

C(t w + T,t w ) = ^2q j [U j (t w + T,t w )-Il j+1 (t w + T,t w )] =^2[qj-qj-i]Uj(t w +T,t v) ) , (36) 

j=0 i=0 

where j labels the level of the tree from top to bottom, and Hj(t w + r, t w ) is the probability 

that no jump beyond the j level of the tree has occured between t w and t w + r. All the 

levels j > M* such that x(qj) is larger than 1 are equilibrated on microscopic (~ t ) times. 

This means that the corresponding Hj are zero as soon as r > t„. This part of the tree thus 

contributes to the stationary dynamics, while the levels j < M* contribute to the aging 

part 

M* 

C(t w + r, t w ) = Y^ Qj[Kj(tw + r,t w ) - Uj+i(tu> + r, t w )] + C ST (r) . (37) 

j=o 

It is easy to show 113 that the short time (r <C t w ) decay of the aging part of C(t w + r,t w ) 
behaves as in Eq. (30), with x — ► x(qM*), while the long time (r 3> t w ) decay is decribed 



■'With the convention ITm+i = 0, and q_i = 0. 
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by Eq. (pTI) with x — ► a; (go)- Not surprisingly, the short time decay is mostly sensitive to 
the fastest part of the tree j = M* , while the long time decay is governed by the slowest 
processes j = 0. Note that the experimental data is indeed such that the long time x(qo) is 
smaller than the short time exponent x{qM*)- 

The presence of levels such that x ~ 1 is very interesting from the point of view of 
'l//'-noise, for the following reason. If x > 1, the corresponding contribution to the a.c. 
susceptibility is stationary and behaves as XstO 4 ') ^ ( UJ t ) x ~ l . The noise spectrum is thus 
given by 

2T X "M 2 

b{LU) = oc uj . {6a) 

IT UJ 

So the leading contribution to the low frequency noise coming from these equilibrated 
processes with x > 1 comes from x — 1 and scales as 1/u. On the other hand, if x < 1, the 
contribution to the a.c. susceptibility is aging and behaves as \"{u),t w ) oc (wtu,)^ -1 . This 
contribution is thus decaying with time, but more and more slowly as x approaches 1 from 
below. Among these aging modes, those which are the slowest to disappear correspond 
to x = 1", leading again to a 1/u dependance. Within this picture, 1// noise appears 
rather naturally; furthermore one expects that this noise should generically exhibit some 
non-stationary contributions. 

Another interesting aspect of the dynam ics o n a 'multi-level' tree is its response to 
temperature cycling. It was suggested in Refs[ 13 IPl] that the negative temperature cyclings 
- which reveal both 'rejuvenation' in the intermediate, low temperature period, followed by 
a perfect 'memory' of the a.c. susceptibility - point towards a hierarchical picture of phase 
space, where finer details are progressively revealed as the temperature is lowered. Using 
the fact that the whole curve x(q) decreases when the temperature is lowered, it is easy to 
account for this phenomenon within the multi- level tree model E3. 

• Energy or 'entropy' barriers ? 

Before discussing the possible relation between this hierarchical picture and real-space, 
droplet like descriptions, one should emphasize that in the trap models discussed above, 
aging is induced by the presence of energy barriers, the crossing of which becomes slower 
and slower as the temperature is decreased. Aging is nevertheless also present in models 



where there is no barrier crossing at all. It was already noticed 1123 that the SK model at 
T = has a slow decrease in energy and never reaches a stable configuration. This same 
situation can be seen more clearly in the following model: let us consider the casea where 
the hopping rate W a ^p is equal to zero if the state energy fp is larger than f a , and equal to 
Wo = (Nto)^ 1 otherwise. This rule actually corresponds to the Glauber dynamics at zero 
temperature. In the limit N — ► oo, the system never reaches the ground state: there always 
exists states of lower energy towards which the system can evolve - the number of these 
'escape directions' however becomes smaller and smaller as time increases. One can show 
that independently of the distribution of energies p(f) (but provided that these energies are 
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independent), the correlation function in this model is given byE3 

C(t w + r,t w ) = (q EA -qo)— -^ h<?o- (39) 

t w +T 

Note that the decay of C(t w + r, t w ) is regular when r — ► 0, in contrast with the above 
trap model and, as we shall see, with the generic mean-field situation. Slow dynamics in 
this model can thus be attributed to 'entropic barriers', i.e. the fact that paths leading 
to smaller energies become more and more scarce as time increases. A similar scenario 



holds in the 'Backgammon model' introduced by Ritort 11211 and studied in detail in Ref. 
p23]. It is reasonable to expect that both energy and entropy barriers should contribute to 
non-equilibrium dynamics in real systems. 

• Speculations about trees and clusters. 

What could be the interpretation of Eq. ( p6| ) in real space ? Clearly, large q's should be 
related to small clusters of reversed spins, corresponding to 'fast' processes, while small q's 
correspond to large clusters, or 'slow' processes. Let us suppose that the overlap between 
two configurations can be written as a sum of contributions from 'clusters' of different linear 
scales £: 

1 L i^Y 

1*3 = yiEE QtQijjih) , (40) 

1=1 u=l 

where L is the size of the sample, qe is the incremental contribution to q of the clusters of 
size £; for fractal clusters]] of dimension dj one expects qi ex £ f . ® a n(it) is equal to zero 
if the states a and f3 differ by the reversal of the cluster of size £ in the 'cell' labeled it and 
equal to one otherwise. Following the same speculative vein, one can write the two-time 
correlation function as 

L 

C(t w + r,t w )=J2% u e(tw + T,t w ) (41) 

1=1 l 

where He(t w + T,t w ) is the average (over all the 'cells' of size £) fraction of clusters which 
have not flipped between t w and t w + r. The disorder average is no longer needed here since 
there is a space average over many independent clusters. 

Assuming that the barrier heights are distributed exponentially with a scale depen- 
dent parameter x(£), we find that the correlation function will behave much as above in 
Eqs.(30,[3l]). The parameter^!/) actually fixes the relation between energy scales T/x(£) 
and length scales (see Refs.pS]). Taking x(£) = T/(T£ e ) (as su gge sted by scaling argu- 



ments or replica calculations on the problem of pinned manifolds l 104 lr 9 lr "1) , one finds that 
there exists a characteristic length scale £* such that x{£*) = 1, separating small length 
scales £ < £* - for which equilibrium is reached - from large length scales where aging takes 



fc The idea of fractal clusters of spins in spin-glasses dates back to Ref I23i\. 
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place. In other words, 

L 

C(t w +T,t w ) = J2% U e( t w+T,t w ) + C st (t) (42) 

£=£* 

with C st (t) = J2e=i ffn^(r). m particular, one has: 

L 

q EA = lim lim C(t w + r,t w ) = V -7 . (43) 

T->oot„^oo tte*t 

Assuming thatH T oc (T g - T) w , one thus finds q EA oc (T g - Tf, with = ((d- d f )uj)/6. 
Interestingly, as the temperature is decreased, there is an infinite sequence of 'glass tran- 
sitions', where all the length scales (in decreasing order) are progressively driven out of 
equilibrium. 

• Conclusion 

Although there is still a lot of work to do to clarify the above picture and make it 
consistent, the idea of modelling the dynamics of a complex system through the motion 
of a point 'particle' in a random potential is fruitful, and actually used in many different 
contexts (structural glasses, protein folding, etcH). The aim of the present Section was to 
show that such models can naturally lead to aging. Actually, the next Sections, devoted 
to an analytical study of some mean-field models of spin glasses (which can also be seen 
as models of diffusion in a random potential) will share many similarities and important 
differences with the above discussion. 

3 Mean-field models of aging: analytical results 

It took several years to realize that mean-field models of spin- glasses, endowed with a 
suitable relaxational dynamics (usually Langevin, though Glauber is also possible), actually 
do capture some aging phenomena in the glassy phase cSeII'Ej. They thus provide a set of 
microscopical models where glassy dynamics and aging effects can be studied analytically. 
This Section will summarize the main results obtained in the recent years on these mean- 
field models, which have been most valuable in clarifying some of the basic theoretical issues 
described in the Introduction. The relation between these models and finite dimensional 
systems is still very much a matter of debate; we shall however postpone this discussion to 
Section 3.8 and the Conclusion. 

The basic simplification occuring in mean-field models is that, after averaging over the 
disorder and making the number of spins large (N — > 00), one obtains a set of closed 
equations for the two-time correlation and response functions, from which the energy and 
magnetization can also be calculated. As we discuss below, these equations imply the exis- 
tence of a critical temperature T c , below which aging effects appear, and the fdt is violated. 
On the other hand above T c , the same equations allow for a tti solution, compatible with 

FDT. 
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Most mean-field dynamics studied so far have focused on models which belong to the 
following family of spin glass Hamiltonians, describing the interactions of N continuous 
spins (j>i, i = 1...N: 

oo 

E(W})=J2 F r E Jh,i2,...,ir+l<t>il-<f>ir+l (44) 

r—1 ii<i'2,...<i r +i 

where the J,-, ,■„ ,- ., are random Gaussian variables with variance A r ~' r ~ 1 . Different choices 
of F r lead to different models. 

A quartic spin weight term can be added in order to make a 'soft' version of Ising spins, 
or one can consider a spherical version, i.e. J^ (f>? = N. A popular choice is the p-spin 



spherical modeled defined by F r = g5 r+ i )P , plus a spherical constraint. In what follows we 
shall mainly use as an example the spherical (or Gaussian) versions. The same model (Ell) 
can also be seen as the potential energy of a point particle in a random potential, where the 
4>iS are the coordinate of the particle's position in a N dimensional space E£j. In this case, 
the energy is a Gaussian random potential, the correlations of which are related to the F r 's 
as 

E{{4>})EM}) = NV f - £ M V(x) = £ j^fjy* r+1 (45) 

where the over line means an average over the quenched disorder. 

We shall in the following consider the dynamics to be modeled by a Langevin equation: 

^ = -iA(jt)<f>i -qT: + Vi(t) + hi(t) (46) 

where the white noises rji are mutually independent and of variance 2T, and hi(t) is a 
time dependent external field. The 'mass' term n(t) is incorporated in order to enforce the 
spherical constraint each time, or to model the presence of an harmonic potential in the 
case of a particle in a random potential, but may be set to zero in other cases. 
The correlations and response are defined as: 



, N -. N 



s(M*)) 



^x^w-r.v u, V , N Z^ §h . {f) 



(47) 



h t =0 



where the braces mean average over the thermal noises r\. There exist by now well estab- 
lished methods in order to obtain the equations of motion for these systems in the large 
N limit (time being kept finite). The best known is to introduce a dynamical field theory 
partition function, average over disorder the partition function by using the fact that it is 
normalised 1153, and compute it for large N by a saddle point method. This is the route 
which was followed originally by Sompolinsky and Zippelius&tl. An alternative derivation 
uses thaxavity method. We refer the reader to the original papers or to more recent text- 
books qH, and rather focus on the solution of these equations. Starting the dynamics at 
time t = from a random configuration (chosen with a uniform distribution in configuration 
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space, corresponding to a quench from infinite temperatures) the dynamical equations for 
the spherical or the Gaussian case are found to be 

^M = -n{t)C(t,t')+2TR{t',t) 

f rt 



+ f ds D(t, s) R(t', s)+ [ ds E(i, s) C(s, t') , (48) 

JO Jo 

-Hit) R(t, t') + Sit -t')+ ( ds E(t, s) R(s, if) , (49) 

Jt' 



/o Jo 

dR(t,t>) rt 

di 



where 



E(M') = R(t,t')V"[C(t,t')] , 

D(t,t') = V'[C(t,t')}. (50) 

It is worth keeping in mind the limitations of the present approach. For example, the simple 
form for T,,D in terms of C, R is a peculiarity of this class of models. More complicated 
forms in which £(£, t'),D(t, t') are Junctionals of C, R (with integrals involving C and R at 
intermediate times) rather than ordinary functions are obtained, for example, in the case of 
the Sherrington-Kirkpatrick model. So long as the functional dependency is only on C, R, 
one may however expect that they can be treated with the same methods. On the other 
hand, as soon as one introduces a Hamiltonian with a finite number of neighbours per spin 
(even for mean-field-like Hamiltonian such as Bethe lattice, or random lattice systems), the 
dynamical equations do not close on the two point correlation and response. One must 
then introduce a whole hierarchy of k point correlations and responses, which has not been 
investigated yet. 

3.1 Self-averageness and 'universality' 

Equations (|48|) , (fl|) and (|50|) are exact in the large ./V limit — for times that do not diverge 
with N. Furthermore, one can show using the same methods by which they are derived, 
that the correlations and responses are self-averaging with respect to both the thermal noise 
and the realization of disorder, again for times that do not diverge with N. Hence, we could 
well have omitted the braces and the overbar in Eq. ([17]). Non self-averageness in certain 
macroscopic quantities appears only for times that diverge with the system size, and in 
particular within the equilibrium Gibbs-Boltzmann measure. 

The insensitivity of dynamics with respect to the realisation of disorder is probably 
intimately related to the fact, discussed in the next Section, that certain models without 
quenched disorder show very similar out of equilibrium behaviour, and hence can be studied 
by considering them as disordered. The correspondence between the random and non- 
random version might however break down for (divergent) times when non-self averaging 
features appear. 



24 



8 



0.(3 



0.4 



0.2 



1 










\ ' 


1 


1 


0.8 










0.6 










J 








- 


0.4 










0.2 








FDT 



t 


, 


, 


, ^\ 


\ 


0.2 


0.4 0.6 


0.8 










y 



present 



Figure 6: The response R(t, t') versus t' — t at high temperatures T — 1 > T c for four total times t = 

ti,t2,ti±jt4. The curves were obtained from a numerical resolution of the dynamical equations (|4q), (W9) 

and (BOf) in a p = 3 spherical spin glass. All four curves merge into one, R(t,t') = R(t — t'). In the inset, 

X(t,t') = x(t — t') vs. C(t,t') = C(t — £'), fdt is satisfied and the slope is just —1/T. 



3.2 The high temperature phase: two types of spin glasses 

One expects that the above spin glass models converge fast towards a paramagnetic equilib- 
rium phase at high enough temperatures, where the system has no long-term memory and 
obeys tti and fdt. This can be seen nicely from a numerical study of the dynamical equa- 
tions at relatively short times, which will also be useful to identify the qualitative behaviour 
at low temperatures. In order to study the memory of the system, we plot the response 
functions R(t,t') versus t' — t at different values of time t, t = £i,£2>*3>*4- (These curves 
were obtained numerically from the model model defined in Eq.(|44|) with only F2 7^ 0. A 
simple numerical procedure consists in discretizing time evenly and iterating the dynamical 
equations which are causal. With some extrapolation procedure on the mesh of the grid 
one can reach safely times of order 1000 Ej.) 

Figure || shows that the response at high temperatures does not extend to the distant 
past and that it depends only on time differences. The system achieves equilibrium after 
a transient, and eventually forgets the origin of times. Actually, the correlation function is 
also tti. Finally, the x vs. C plot (see Section L3) is a straight line of slope —1/T for all 
values of C (see the inset in Fig. ||) showing that fdt is satisfied. 

Turning now to the analytical treatment of the equation, it is easy to show that if tti 
holds (i.e. if the two unknown functions of two variables C(t w + r, t w ) and R(t w + r, t w ) 
actually depend only on r), then the second equation ( |49| ) is a consequence of the first ( [48| ) 
provided the fdt is satisfied. We are thus left with a single function C(t) in Eq. ((481), and 
the kernels £ and D are related by dD/dr = — S/T. Using the finite extent of the memory 
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we can safely send the initial time to — oo, to obtain: 

^ = -^ooC(r) + ^f_Jr'V[C{r - r>)]^± , C{r = 0) = 1 (51) 

where 

Hoc = K* -» °°) - V[l]/T . (52) 

The above equation is valid as long as C(t) decays to zero in the long r limit. 

Equation (|5l| ) is basically the general 'schematic' Mode-Coupling equation for the den- 
sity correlations in a supercooled liquid above the dynamical transition temperature intro- 
duced by Leutheusserli^l, Gotze and others Ii3 as a model for the ideal glass transition. The 
only difference lies in the fact that the Mode-Coupling equations also possess an 'inertial' 



term d^.C(r). This coincidence will be further discussed in Section 4.4 . 

The behaviour of the solution to these equations when one lowers the temperature 
depends upon the structure of the disorder, i.e. upon the function V. It turns out that 
there are two broad classes of mean-field spin glasses, characterized by rather different 
behaviours. The point is that as one lowers the temperature, there appears a critical point 
T c at which the decay behaviour of C(t) shows a marked change. For 'discontinuous' 
models, C(t) does no longer decay to zero, but rather to a finite value q EA > (which 
is also called the non ergodicity parameter / in the context of glasses). For 'continuous' 
models, C(t) still decays to q EA = at T c but with critical slowing down. Below T c , q EA 
grows continuously from zero j. 

• 'Discontinuous' models 

The simplest prototype, which we discuss here, is the pure spherical p-spin model with 
F r = <?cV+i,p- Another example is provided by the problem of a particle in a random 
potential with short range correlati ons i.e. when V{x) decays sufficiently fast for large x's. 



The solution, studied in Refs. f^ l 108 ^ 09 ! ] behaves as follows. Above the critical temper 
ature T c , the correlation C{t) decays to zero at large r. Slightly above T c , the correlation 
already starts developing a plateau at C ~ q EA before eventually decaying to zero, as shown 
in Fig. 0. The length of the plateau increases as a power law of T — T c when the temper- 
ature gets closer to T c . The details of how C{t) first decays towards the plateau and then 
departs from it has been one of the most studied aspect of the mct in the context of glass 
forming liquids, since these features can be directly tested experimentally using a variety 
of techniques. One finds that: 

C(t) ~ q EA + c a T~ a C>q EA 

C(t) ~ q EA -c b T b C<q EA 

where the exponents a, b are related by 

r 2 [l+6] _ r 2 [l-a] _T c V"{q EA ) 

F[l + 26] T[l - 2a] 2 (V"(g EA )) 3 / 2 ' l ' 

'These two types are also called, respectively, A and B in the context of the Mode Coupling theory E3. 
One should keep in mind that we speak here of dynamical phase transitions, and this terminology is not 
related to the Ehrenfest classification of equilibrium phase transitions. 
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Figure 7: The correlation for the same model as in Fig. q at high temperatures close to the transition T ~ T c 



Clearly, the equilibration time for a system that is quenched to T > T c is at least as large as 
the time of the plateau t(T). A formulation in terms of time-differences ( |52[ ) is thus valid 
only when the waiting time is much larger than the a-relaxation time t(T). For times that 
are smaller than t(T), one has to go back to the two-time equations — this will be the case 
throughout the low-temperature phase. 



• 'Continuous' models 

The class of 'continuous' spin glasses contains the more usual case of the SK model El 
(which is however not described by an equation of the type (^)). Among the systems we 
are discussing is the case of a particle in a long-range correlated random potential t!3 (V(x) 
decaying as a power law), or some spin systems such as a mixture between p = 2 and p = 1 
interaction 



Hi. 



F r — g6 r +i 2 + 9 <Wi,4- 



Again there is a critical temperature T c above which the correlation C{t) decays to 
zero at large r. The main difference with the previous case is the absence of the plateau 
structure around q EA for T slightly above T c , which is obviously related to the fact that in 
this case the Edwards- Anderson order parameter departs continuously from zero when one 
decreases the temperature through T c . 



3.3 Low temperatures: Weak long term memory and weak ergodicity breaking 

Let us now discuss what happens below T c . Again, we first show the numerical solution 
of the full causal dynamical equations ([48|), (|49l) and (|50|), and look at the same plot as 
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Figure 8: The response R(ti,t') in terms of t' — t% at a low temperature T < T c for the same model as in 
Fig.M. For each curve t, is t\ = 100, ti — 200, £3 — 300, £4 = 400, respectively. 'Quench' corresponds to 

t' = while 'present' to t' = U. 

before, namely the response functions R(t, t') versus t' — t, for different values of time t, 
t = ti, t2, £3, £4- Figure || shows that: 

- The system has a strong response to pertubations in the immediate past that is quite 
similar to the high temperature response. 

- However, a long tail extending down to the quench time t' = has now appeared. 
The total area under the response curves, f ds R(t, s) approaches at large times a finite 

limit, which is equal to the linear susceptibility x to a constant field. Part of this area is 
already given by the peak to the right of Fig. S, which is the high-frequency, stationary 
contribution to the susceptibility. It turns out however that the area below the long-time 
tail also gives a non-zero contribution - the memory to the distant past is substantial and 
can never be neglected: 



Vt* 



lim / dsR(t, s) < lim / dsR(t, s) 

:^oo J t _ t * t^oo J 



(55) 



This can be hinted from the simulations. As we shall see, it can also be derived from the 
dynamical equations. 

At this stage it would seem that in order to solve for the correlations and responses at 
large times we need to know the complete solution at all times, because the memory kernels 
(which involve the response) have nontrivial contributions from all the past. If this were 
the case, the problem of finding an asymptotic solution for long times would be hopeless 
from the analytical point of view! 



28 



c 




quench 



present 



Figure 9: The correlation function C(ti,t') vs t' — ti at low temperatures for the same model as in Figs |6| 

and |. From left to right U = 400, 200, 100, 50. 

Fortunately, something remarkable happens within the models considered here: even 
though the area under the long-time tails of the response remains finite, the height of the 
tails themselves tend to zero as we consider larger times. More precisely, the integrated 
response at time t to a signal between the initial time s = and any finite time s = t* tends 
to zero at large t: 



lim 

t— >oo 



dsR(t,s) = 



(56) 



for any fixed t*. The same can be said about the integral of R x G where G is any finite 
function, such as the memory kernel S. This behaviour was described as 'weak long term 
memory' in Ref. p]: the memory tends to be 'weak' for any finite interval of the past, 
but is strong when integrated over the whole past FT. This is close in spirit to the 'weak 
ergodicity breaking' property defined in Section ( |1.4|) : a perturbation lasting for any finite 
duration will eventually be forgotten. Therefore the long time dynamics will decouple from 
the initial (non universal, and out of control) transients. 

The evolution of the two-time correlation function also has (at least) two distinct 
regimes, as shown in Figs. |9| and [[(]. For times such that r = t — t' is small all curves 
merge and one has tti. (These are times close to 'present' in Fig. |9| and to the left in Fig. 
|l0| .) The corresponding decay is 'fast' (see footnote 6 above) in the regime where C drops 
from 1 at equal times to the plateau value q EA , defined in Eq.(|]). However, when C decays 

m Not every system will satisfy this condition; on the contrary, systems that remember their initial 
transients are much harder to treat and how to deal with their dynamics is still a fully open problem. 
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below q EA it does it in a manner which depends on both t and r. This subsequent decay is 
'slow'. 

This suggests (and the analytical calculation later confirms) that one can perform, for 
both the correlation and response functions, a decomposition into a stationary and an aging 
part, similar to the one introduced above for the description of aging experiments (see Figs. 
|and||): 

R(t w + r, t w ) = R AG (t w + t, t w ) + .R S t(t) , 

C(t w + r, t w ) = C AG {t w + r, t w ) + C st {t) . (57) 

Defining as in Eq.(g) the Edwards-Anderson order parameter: q EA = linv-,.00 limt TO -^oo 
C(t w + T,t w ), the stationary (and thus the aging) parts are defined by: 

C st (t) = lim C(t w + r,t w ) -q EA R st (t) = lim R(t w + T, t w ) (58) 

t w — >00 t w — >OQ 

It turns out that the stationary parts, which are by definition tti, also satisfy the fdt, 
Rst( t ) = —j;dC ST (T)/dT. Finally, note that the above definitions imply: 

lim^^oo C AG (t w +T,t w ) = q EA lim^^oo R AG (t w + T,t w ) = 

linx™ C st (t) = C ST (0) = l-q EA (59) 

linv^oe i? ST (r) = R ST {0) = 1 



3.4 Low temperature solution of the dynamical equations 



• General strategy. 

An asymptotic solution to the dynamical equations was first fou nd in Ref. pil] in the 
p-spin spherical model and then gener alise d to various situations c3~Eitij'tiil. It has been 



described in detail in several papers laOEil. We restrict here to the essential assumptions 
and the main ideas allowing one to find the low temperature solution of the dynamical 
equations @j and (H). 

The solution strongly relies on the weak long term memory assumption. It is asymptotic 
in the sense that it holds only at large times t w : the transient effects are much more 
complicated, and have not been studied yet. The method of solution we shall outline here 
only determines the time-dependences in the aging regime up to a reparametrization in time 
t — > h(t). Thus, one obtains a family of solutions for the aging part, related to one another 
by time reparametrisations: 

C AG (t,t') - C AG (h(t),h(t')) R AG (t,t') -> (^P) R AG (h(t),h(t')) . (60) 

(Note that the presence of the factor dh(t')/dt' comes from the fact that it is the integral of 
R over time, rather than R itself, which is the physical quantity). The 'selection problem' 
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of determining the actual function h chosen by the system is still an open one that requires 
the matching of the regimes of short and long time differences. 

The starting step consists in proposing asymptotic forms for the aging part of the 
correlation function in the two-time plane valid for large t w + r and r. As we mentioned in 
(|1.2|), a possible form for the correlation is: 



c AG (t w + r, t w ) ~ y]c 



i 

i 



tli\t b 



hi(t w + t) 



(61) 



where each term will vary in separate time sectors, defined by taking the times to infinity 
with < hi(t w )/hi(t w + t) < 1 (see Eq.([lO|)). The asymptotic form of the correlation is 
given then by the knowledge of Ci and hi . This form is meant to represent the correlation 
only in the limit of large times, and it need not be unique. One may wonder whether it 
exhausts all possibilities. 

This question can be answered with the following construction^. Consider the con- 
figurations at three large times t m i n < t; mt < t max , and the corresponding correlations 
C(t max ,t min ), C(t max ,t int ) and C(t in t,t min ). Using the fact that the correlations decrease 
with time-separations, one can show that in the limit of large times the three correlations 
must be related by 

^ J [tmax,tmin) — J \S^ [tmax, tint), ^ (tint, tmin)} v" / 

where the function / defines the geometry of the triangles described by the trajectory of the 
system in phase space. Now, it is easy to see that / is an associative function, and one can 
classify all the possible forms of an associative, monotonical functioned using elementary 
group theory. 

This in turn leads to a classification of all the possible two-time scalings as follows. One 
considers the special (fixed point) values of the correlation qo, q\, . . . , qk = qEA, Qk+i = 
C(t,t) defined by f(qi,qi) = q%- If the number of fixed points is finite, one can prove 
that the correlation can be written for large times as (|6l| ) with qi = J2j<i^j{^)- The in- 
termediate values of correlation between two fixed points constitute a 'correlation scale', 
within which only one term of ( |61~D is non-constant. Triangles whose sides belong to dif- 
ferent scales (e.g. C (t max , t int ) < qi , C(t in t,t min ) > %) are isosceles with C(t max ,t min ) = 
uiin[C (t max , ti n t), C (ti n t, t m i n )]: there is ultrametricity tiJ between correlation scales. It 
may also happen that the fixed points qi form a continuum within a certain range (q a , qb)- 
In that case we have, for any two correlations within this range: 

(- / \tmax,tmin) — mm|G [tmax, tint), ^ [tint, L"inin)\ ^S t m in * OO . [P^l 

Correspondingly, the correlation has to be represented in that case as a limit of a continuous 
sum of infinitely many scaling functions Ci with vanishing weight. Note that the above 
construction does not rely on any concrete model and is not restricted to mean-field. 

It is interesting to see that this form of ultrametricity in the correlations appears 
in a very natural way within out of equilibrium dynamics. It is already present in the 
simplest form in the case of domain growth: if we consider three very large times such 
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that C(t max ,ti n t) > rn 2 and C(ti n t,t m i n ) < m 2 (where m is the magnetization), then 
C(t max ,t min ) = mm[C(t max , ti n t), C (ti n t, t min )]. Ultrametricity between correlations larger 
and smaller than m 2 expresses the fact that the time scale corresponding to relaxation due 
to domain wall motion becomes infinitely larger than the one corresponding to thermal 
fluctuations within each domain. It is however not clear that there are finite dimensional 
models for which a 'full' ultrametric form (|63| ) holds for any two values of the correlations. 
Having discussed the possible asymptotic forms of the correlations, the next step is to 
make an Ansatz for the response function. One can write without loss of generality: 

X(t,t')dC(t,t') 
R(t,t) = — ^— (64) 

(t > t'). The Ansatz now consists in proposing that for large t,t', X only depends on time 
through C: 

X(t,t')=X[C(t,t')] (65) 



or, in other words, that the parametric plots x vs C in Section |1.3| converge to a limit 
curve as t w — > oo. The asymptotic form for the response is thus obtained from that of 
the correlation through the introduction of a certain function X[C]. At this point one 
substitutes the above Ansatze for the correlation and response functions into the dynamical 
equations. In this way one determines 

i) X[C] (which contains the information on aging of the response) and 

ii) /, or, equivalently, the number of terms in (pll) and the respective Cj. 

Hi) One should in principle also obtain the hi, but this requires to solve the selection 
problem discussed above. 

The corresponding computations can be rather lengthy. We shall not detail them here, 
but rather present the results obtained when one applies this technique to the mean-field 
models considered above. 

• The case of discontinuous models. 



For low temperatures T < T c the correlation is given by Eq.(61) with only one aging 
sector (i.e. one function h(t)) plus the tti scale. It reads 

C(t w + r, t w ) = C st (t) + C AG f ,^ ( *7 ) r ) ■ (66) 

\h(t w + t)J 

The memory properties of the system are controlled by X[C], which takes a particularly 
simple form: 

X[C] = 1 for C > q EA (fdt) 

X[C] = X < 1 for C < q EA . (67) 

A is a positive, temperature-dependent number, constant throughout the aging regime. 
The fact that X is constant is not assumed a priori, but comes out from the equations of 
motion. A possible explanation for this fact is found when interpreting T EFF = T/X as 
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an 'aging temperature^: it means that degrees of freedom of comparable frequencies are 
mutually thermalised e3. Note that at the glass transition, the aging temperature is equal 
to the bath temperature. The ratio T EFF /T increases with decreasing bath temperature; in 
particular, T EFF remains non-zero at zero bath temperature. 

The behaviour of the correlation around the plateau q EA provides the low-temperature 
extension of the one already encountered at high temperatures. Taking the limit of large 
t w before taking the limit of large r one only explores the approach of the correlation to 
the plateau q EA . This decay is still given by a power law characterised by a temperature- 
dependent exponent a, as was the caseJor T > T c . However, the subsequent departure from 
q EA is t^-dependent and characterized e3 by another temperature-dependent exponent b: 

C(t w + T,t w ) ~ q EA + c a T~ a C>q EA 

( \ h < ( 68 ) 

C(t w + T,t w ) ~ q EA -c h [^j C^q EA 

where T w is an effective waiting time, defined as T w = h(t w ) / h! \t w ) . (In the case in which 
h{t) is a simple power law, one has T w oct w .) The exponents a,b have precisely the same 
meaning as the exponents a, b defined in Section ([hi]), and are in this case related by 

Y r 2 [l + 6] r 2 [l-q] = T V"(q EA ) 
r[l + 26] r[l-2o] 2 (V"(g EA )) 3 / 2 l ' 

with q EA given by Eq.(J6|). 

• The case of continuous models. 

As regards the temporal behaviour of the correlations in the stationary regime, the 
approach to the plateau at q EA is also given by a power law with a temperature-dependent 
exponent a. The situation within the aging regime is however more complicated than in the 
case of discontinuous models. The aging part of the correlations satisfy the full ultrametric 
triangle relations (|63|): 

f(C 1 ,C 2 )=mm[C 1 ,C 2 ] (70) 



if at least one of C±, C 2 is smaller than q EA . A representation like (p8[) for the escape from 
the plateau is not possible for these models unless one makes the exponent b waiting-time 
dependent. 

For purely continuous models the function X[C] is not a constant and, remarkably, 
coincides with the function x(q) in the replica treatment of equilibriumj]. This and other 
coincidences between out of equilibrium dynamics and statics of continuous models have 
escaped any kind of physical understanding so far. (For some recent work in this direction, 
see Ref. 



"The dynamical X[C] is much easier to obtain numerically than its static counterpart because one does 
it have to equilib 
astonishingly goodt 
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not have to equilib rate, i The numerical confirmation of the analytical form of X[C] for the SK model is 



3. 5 Generalization to several coupled modes — the case of spatial dependence 

New physical insights appear when we consider the generalization to several coupled models, 
and, in particular, to mean-field cases in which there is also spatial dependence. In such 
cases we have to deal with several correlations Ckk'(t,t') and responses Rkk'(t,tf), where 
the indices k, k' refer to different modes, for example they can represent spatial positions 
Xk or Fourier components. 



The construction of an Ansatz proceeds as before 112 !. One has to add now a prescrip- 
tion for the long-times relationship between different correlations. Choosing one particular 
correlation function Co.o as an effective 'clock', one may look for solutions of the form: 

Ck,k>{t,t ) = Fk,k'[Co,o{t,t )] (71) 

with Tk,k' to be determined. One also introduces the fluctuation-dissipation violation factors 
defined by 

RwW) = x wl c wM)}Ww . (72) 

Interestingly, it turns out that the Ansatz closes with two extreme possibilities: 

i) Xkk' 7^ f° r k ¥" &'i an d ^k,k = Xy y at equal times; 

ii) Xk^k 1 — ► for k / k', and Xk,k, Xy t y possibly different. 

This can be understood as a property of partial thermalisation c3: remembering that 
Tk = T/Xki- is an effective temperature, in the case i) the subsystems have T^ = Ty at 
corresponding time scales, while in case ii) the subsystems have zero cross response and T\~ 
may be different from Ty. 

A mean-field case with many modes is obtained when one studies the dyna mics of a 
random manifold in a disordered medium within the Hartree approximation 1113. One has 
all the time scalings described so far for each mode k, plus dynamical scalings in terms of 
k and times. It turns out that a solution satisfying (fnj) and (|7^) appears naturally in that 
case. 

3.6 Speculations on the 'effective' age function h(t) 

Although the above solution describes some general features of the low temperature, aging 
regime of mean-field spin glass models, it is still, even in the simpler discontinuous case, 
incomplete. As mentioned above, the dynamical equations become, in the asymptotic (t — > 
oo) limit, invariant under any monotonous time reparametrisation t — > h(t). The function 
h(t) can in principle only be determined through matching with the early (r ~ t ) solution, 
or by a numerical solution of the two-time equations. For the spherical p-spin problem, the 
latter procedure suggests that h(t) is close to a power- law cH, in other words that the aging 
part of the correlation function C(t w + r, t w ) is a function of the ratio t jt w . In principle, 
h(t) could be any other function of time, for example: 



h(t) = exp 



1 /tX 1 ^ 



1 — M Vt 



or h(t) = exp 
31 



log^(f) 



(73) 
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in which case, the effective 'age' T w appearing in, e.g. (|6cf ) is 
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(74) 



Except in the cases /i = 1 or v = 1, where /i(t) = t/t , these more complicated forms lead 
to an explicit dependence of the effective age T w on the microscopic time scale t . (More 
precisely, when \x < 1 (or v > 1), the effective age T w is much smaller than t w , a situation 
called 'subaging' in Ref.p].) 

The simplest scaling situation would be that the value of t becomes irrelevant on the 
experimental time scale where t w ,T ^> t . This would lead immediately to 'full' aging, 
T w oc t w . Naive scaling can however break down in some cases ti3, which means that the 
value of t is important even in the t w — ► oo limit - in other words that the effective age 
T w depends on t . Such is the case of systems with logarithmic domain growth e.g. the 
random field Ising models, where T w oc t w \og{^-) ('superaging'). Signs that subaging may 

We feel that the 
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also happen in the mean-field models considered here can be found in 
determination of h is one major unsolved issue in mean-field dynamics. Correspondingly, 
the same ambiguity remains from an experimental point of view: a detailed analysis of the 
TRM reveals small but systematic deviations from 'full' aging. A scaling function hit) of 
the above form, with [i =jQ.97 or v ~ 2, does a better job at fitting the experimentslij (see 
Figs. 3.b and 3.c. in Ref.p]). However, these deviations might alternatively be interpreted 
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being 'interrupted aging' (i.e. equilibration in a long but finite time) due to finite field 



1C ,118 



or finite sizeSLIj effects. 
At any rate, note that a pure r/t w behaviour would rule out the existence of many time 
sectors with full ultrametricity, and would leave a scaling as in the discontinuous models as 
the only possibility. 

3.7 Out of equilibrium versus 'equilibrium on diverging time scales' 

Dynamical studies of mean-field spin glass models in their low temperature phase are not 
w. They started more than 15 years ago with the work of Sompolinsky and Zippelius 
I and were subsequently applied to many other problems. The above description, which 
focuses onto out of equilibrium dynamics, has followed a very different route - which was 
only found recently ei. However, as all these works address similar issues, a comment on 
their relationship is in order. For lack of space we shall not be able to present in any detail 
the former approaches, but these are by now well documented EkI'BqQ. We rather want to 
stress the conceptual differences with the present approach. 

Consider a mean-field spin glass model below T c . The two times dynamical equations 
(|4q) and ( (49|) are exact equations relating the correlation and response in the mean-field 
models. They describe the behaviour of the system at times that do not diverge with N. 
Even when speaking of 'long times' within this framework, one really means: 

lim lim . (75) 

t,t'-+oo iV->oo 

Obviously for finite N the spin glass will equilibrate in a finite time. There exists an equili- 
bration time t ERG (N), such that for t w much larger than t EnG (N) the system is equilibrated, 
which means that the configurations are sampled with a frequency proportional to their 
Boltzmann-Gibbs weight. The correlation and response then become time translational 
invariant, and related by the fdt. At large N the equilibration time diverges. 

The approach of Ref. p\ starts from the very same dynamical equations (|4q ) and (|49|) , 
but one assumes that the size of the system is finite and large. One also assumes that the 
initial time of the dynamics has been sent to — oo and that tti holds. The corresponding 
construction relies on the hypothesis that there exists a strong hierarchical structure of time 
scales which all diverge with N. Calling t x these time scales, where x is an index chosen 
for instance in [0, 1], the hierarchy means 

lim t x = oo, lim — = oo if x < y . (76) 

N^OO N^OO ty 

The largest of these time scales thus corresponds to the equilibration time £ E rg- 

This allows to produce a solution to the dynamical equations which exhibits a non-trivial 
dynamics within each diverging time scale t x . This dynamics exhibits tti, as assumed 
from the beginning, but it violates the fdt. This astute solution presents several formal 
similarities with the static solution of Parisi (and also with the out of equilibrium study), 
starting with the hierarchical structure, which acquires an appealing interpretation in terms 
of diverging time scales. However it suffers from several problems. At the level of the 
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results first: the correlation function C(t), at very large times r > t ERG , does not go to the 
correct Boltzmann Gibbs equilibrium value, known from the static studies. This is clearly 
inconsistent. A careful study of their derivation reveals two weak points. 

First, there is an inconsistency in the hypotheses. If the initial time has been sent to 
— oo, for a finite N system, then it is fully equilibrated, and then the dynamics is necessarily 
tti and obeys fdt. Otherwise (if the waiting time is much smaller than t ERG ) there is no 
reason to assume either tti or fdt. 

Second, this approach uses the dynamical equations derived within the N — > oo theory, 
for a finite system where activated processes take place. Even though large N saddle-point 
approaches can be used to study activated processes that occur for large but finite N, there 
are subtleties related to the existence of multiple solutions. The multiplicity of solutions is 
related to the fact that the results at times greater than i ERG should be non-self averaging, 
as the replica solution shows ■ 
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Because of these problems, which created many discussions r 2 Tr 2 lr 23 tB, there have been 



various attempts in the litter ature to try to amend this solution, while keeping most of its 



nice mathematical structure. One possibility, first suggested by Horner tin and developed 
in Ref. i23], is to keep the cooling rate finite. This allows to send the volume to infinity 
while keeping a regularization time scale which is the inverse cooling rate, which is sent to 
infinity in the end. What happens once the cooling procedure is over is still not clear in this 
approach. On the other hand, Horner E2J proposed an alternative 'regularizing' procedure 
which consists in making the disorder time-dependent. In this way, aging disappears, the 
solution is tti and manifestly out of equilibrium. 

We believe that the present out of equilibrium dynamical approach, inspired by the 
experiments themselves, is very clear. It is seen to work consistently within the finite time, 
infinite N regime fl75[) . The price to pay is that one has to abandon the postulate of tti, 
and think in the two times plane. In a vague sense, there is also a regularization time which 
is involved, namely the age of the system. 

Finally we want to point out that in all the approaches developed so far to study the 
spin glass phase at the mean-field level, there is a hierarchical (ultrametric) structure which 
is involved. This is true in the statics where it is hidden in Parisi's Ansatz. It is also true 
in the 'dynamics on diverging time scale' approach where it appears in the hypothesis of 
strong time hierarchy (u^). In the out of equilibrium dynamical approach, the situations is 



rather more favourable in that ultrametricity can be proven with mild assumptions (see |3.4 
above). 

3.8 On the links between the static and dynamical approaches. Phase space geometry 

In an infinite system, there is no reason for the out of equilibrium dynamics to be related 
to the static, equilibrium picture. Indeed as we saw, the dynamics refer to finite time scales 
(when N — ► oo), while the static properties are only recovered in the opposite, non-physical 
limit]]. Yet it is instructive to compare the results of the two approaches. It allows to gain 
some intuition on the physical mechanism at the origin of aging. It further underlines the 



°The divergence of tma at large N is known to be of the type £ E rg ~ exp(7V a ). 
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physical difference between the low-temperature behaviour of the two classes of mean-field 
spin- glasses. 

In continuous spin-glasses one gets the natural result that the dynamical transition 
temperature T c coincides with the static critical temperature corresponding to the onset of 
a non zero q EA and of replica symmetry breaking effects. Furthermore, the large time values 
of one-time intensive quantities such as the energy density E(oo) in the limit (|75|) coincide 
with the value E EQ found at equilibrium within the static (replica) approach. 

In discontinuous spin-glasses the results are more surprising: the static transition tem- 
perature T s is lower than the dynamic onaiStH3, T s < T c . In fact the static thermodynamic 
quantities computed within the Boltzmann Gibbs equilibrium are perfectly analytic in the 
neighborhood of the dynamical temperature T c , and conversely nothing special happens at 
T s in the dynamics (|75|). Furthermore, throughout the low-temperature phase, the out of 
equilibrium dynamical energy does not converge to the equilibrium one E EQ . In fact it never 
goes below a 'threshold level' ell E(po) = -Ethres > E EQ . 

• The energy landscape 

A geometric explanation for this strange phenomenon can be found within the frame- 
work of the static mean-field equations of Thouless, Anderson and Palmer (tap equations). 
In the static limit it is possible to write a free energy F TAP (mi, ...,mj\r) in terms of local 
'magnetization' variables which represent the average value of a spin on a large time- window 
(but the large N limit has to be taken first). The minima of this free energy correspond 
to various metastable states; it is known that their weighted sum gives back the correct 
equilibrium results E£3. 

It is far from obvious in general that the dynamical evolution of a system can be seen as 
the relaxation of a point in this free-energy landscape. Yet it is always possible to compute 
the dynamical properties, such as the energy, and to see in what region of the landscape 
the dynamics takes place. This actually provides interesting insights for discontinuous 
spin- glasses. The most complete discussion is available for the case of the spherical p-spin 
system. It turns out that the TAP states can be computed rather easily in this case. The 
reason is the absence of chaoticity: A TAP solution at temperature T\ can be followed 
adiabatically when one changes the temperature to T2, the only change is a global rescaling 
of the (mi, ...,77tj\r). Once the free energies of the solutions are ordered at a temperature 
Ti, this order is maintained at all temperatures: there is no crossing of the solutions in 
the temperature - free energy plane (see Fig. pd| ). Non-trivial tap solutions exist in a 
wide range of temperatures, extending above the static transition temperature T s , and even 
above the dynamic one T c . 

Comparing this landscape to the results of the dynamics, the present understanding of 



the situation is as followaf If H I: above the dynamical temperature T c , there is a coexis- 



tence of a paramagnetic state and a bunch of non trivial, but isolated, TAP states. When the 
dynamics starts from random initial conditions, the system thermalises within the param- 
agnetic state. However if one chooses carefully the initial conditions, one can let the system 
thermalize within one of the tap states, which are separate ergodic components!: 
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Between T s and T c , the paramagnetic state is fractured into exponentially (in N) many 
separate ergodic components, each of which has a higher free-energy compared to that of 
the paramagnetic state, by an amo unt which is exactly equal to their overall configurational 



entropy, called the 'complexity' £9tH-!l. Correspondingly, the simplest static approach (see 
below for more elaborate ones) does not notice this subtlety, and still describes the system as 
a simple paramagnet. In the dynamical approach, starting from random initial conditions, 
one finds that the energy actually remains above the threshold level which is the energy of 
the highest TAP state (see Fig. |ll]). Therefore in the limit (|7q) the system never reaches the 
energies where it would get forever trapped into one of the tap states. This is the origin of 
weak ergodicity breaking in these models. 

In the case of continuous spin-glasses the geometry of the TAP landscape looks superfi- 
cially quite clear: while above the transition temperature T s = T c the free-energy has only 
one minimum, it develops below T c exponentially many states, which matches nicely with 
the Parisi picture. The identification of the low lying tap minima, and of the pure states 
in Parisi's construction to the dynamical ergodic components seems inevitable. The puzzle 
of how to match this equilibrium picture with the out of equilibrium dynamics containing 
infinitely many time sectors is however completely open, as we already saw in the previous 
Section. Surprisingly, many purely static quantities involving different states (which are in 
principle mutually inaccessible) coincide with their dynamical counterparts. In particular, 
the large time values of one-time intensive quantities such as the energy density E(oo) in 
the limit ( |75| ) coincide with the values found at equilibrium within the static approach. This 
means that the set of two dynamical equations ( f48|) and (|4^) contain all the information 
on the replica symmetry breaking solution of these sy stems at low temperatures (A careful 
numerical check of this point can be found in Refs.p'EJ'tU]). Therefore the study of these 



equations may provide an alternative route to a rigorous study of the spin glass phase. 
• The geometrical description of mean-field aging. 

The models we have been describing can be thought of as the motion of a point in 
a high-dimensional phase-space landscape (to the extent that we deal with soft spins). 
For discontinuous models, one finds that, surprisingly, the effect of temperature on the 
dynamics is very minor, provided that one stays within the low temperature phase. Indeed, 
the solution of the equations is regular at T = 0: the dynamics at zero-temperature is 
not very different from the one at finite temperature, T < T c . This at once tells us that 
activated processes are not the main ingredient here. What is then the origin of aging in 
mean-field models? 

At T = 0, we simply have gradient desc ent within the basin of attraction of some 



phase-space minimum, and we can discuss tH3 the essence of the problem without having 



to postulate a 'free-energy landscape' with a dynamical meaning. Since the basins are 
high dimensional objects, a random starting condition will be in the thermodynamical limit 
practically on a border between basins, and hence will remain close to this border for all 
finite times without reaching a minimum. Now, the border is itself partitioned into basins 
of attraction of the stable points on it, which are the saddles with one negative eigenvalue. 
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The same argument as before tells us that, again, they will not reached in finite times: the 
system remains almost on the 'border of a border'. One can now iterate this argument, 
invoking borders within borders and saddles of higher (but smaller than O(N)) indices. 

The conclusion is that high dimensional (phase-space) wells take long to fall into. In 
this sense, the above 'endless descent' scenario concerns equally the motion of a mean- 
field system within a basin of attraction and the 'fall' of a ferromagnet (in any spatial 
dimension) into one of its two 'wells'. Indeed, both or dinary coarsening and mean-field 



aging are reminiscent of the phase-space model described ESI121I in Section 2, in which aging 
is due to the fact that downhill directions are always present, although in decreasing number 
as time grows. What we have described so far is the fact that a macroscopically different 
phase takes an infinite time to grow in the absence of a driving field, a fact that when looked 
upon from the phase-space point of view is recognized to hold even at the mean-field level. 

At finite temperature one cannot invoke a simple gradient descent, but one can still 
use the tap free-energy to understand why the (N = oo) system ages and never falls 
below the threshold level eII. The above discussion can be reformulated as follows: the 
density of eigenvalues of the matrix of second derivatives in each TAP minimum is a shifted 
semicircle law. The smallest eigenvalue A M in(F) decreases as one considers minima with 
higher F, and becomes zero precisely for those with F = F THRES . Above the threshold, the 
spectrum for the saddle points continues to shift, with now A MIN < 0, so that one encounters 
saddles with more and more negative eigenvalues. If one makes a cut of the free energy 
landscape at different values of F, one obtains disconnected 'islands' around each minimum 
for F < -Fthres- As one raises F just above -Fthres one crosses saddles with larger and larger 
number of negative eigenvalues. Each time, a separatrix develops and the set of mutually 
disconnected components becomes more and more connected. One can then picture the 
(N = oo) aging system as falling in more and more disconnected space, hence moving more 
and more slowly - without ever quite stopping, since there are always directions where the 
free-energy decreases (while staying above threshold) ||. 

Again, even at finite temperature this scenario is rather distinct from the 'trap' picture 
of Section 2: Because N is infinite (and the model is fully connected) the system never 
reaches the bottom of a trap from which it could only escape through thermal activation. 
For finite N, however, activated processes will begin to play a crucial role after a finite 
amount of time (see Section ||) . 

• The 'marginal stability' criterion. 

Finally, we would like to mention the possibility of identifying the dynamical transition 
temperature from purely static computations. This is very useful since static calculations 
are generally easier than dynamical ones, in particular when one deals with discrete spin 
variables. 

In continuous spin- glasses, the clearest mathematical characterization of the onset of 
the static transition is to consider two identical copies of the system with the same disorder, 
coupled through a small, but extensive, coupling term of strength proportional to g. For 



p This is similar to the idea of 'percolation' in phase space, investigated in Refs. 

40 



Free energy 
A 




Figure 11: A sketch of the free energy of the TAP states in the spherical p-spin systems. Each TAP state 
like (2) or (4) can be followed adiabatically in temperature until it disappears. The line (1) is the static 
equilibrium free energy, taking into account the multiplicity of the TAP solutions. (5) shows the free energy 
of the TAP states giving the leading contribution to the static equilibrium partition function. The highest 
TAP states (3) are marginally stable and provide the threshold energy: In its dynamical evolution, the energy 
of a system starting from random initial conditions stays above this threshold. 



instance in the Edwards-Anderson model, one can compute the partition function Z 2 for 
two spin systems s and a, coupled through their local energy densities l 130 l: 



H 2 = - 2_, Jij(si.Sj + <Ji<Jj) 



9 ^2 ( Jij s i s j ) ( Jij a i a j ) 

m 



(77) 



The system is in its low temperature phase when the overlap between the two copies (defined 
I / (f3N)d\og{Z2) / dg) is discontinuous at g — ► 0. When g — ► + this overlap tends to 



as 



the Edwards- Anderson order parameter, while when g — » 0~, it is the smallest possible 
overlap between different states. 

In the case of discontinuous spin-glass transitions this construction must be modified 
because of the existence of exponentially many tap states. Roughly speaking the partition 
function can be approximated as 



* = £• 



-l3F a 



/#. 



N(S c (f)-0f) 



(78) 



where a labels various tap states, / = F/N are the free energy densities, and exp(iV S c (f)), 
usually called the complexity, is the number of tap states at a given free energy density. 
The large degeneracy of tap states shifts the dominant region of free energies in the integral 
(|7q) . If one took two copies of the system with a small extensive attraction as in ( |77|) , the 
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leading contribution to Z<i should come from the case where the two spin systems are in 
the same state. However, this leads to Zi{g — * + ) ~ J dfe^ c ^' 2 ^^' which is dominated 
by a wrong saddle point in /. Therefore, in this case, one must be more clever and study 
instead m identical copies of the system in the formal limit m — > 1 (in which case, the 
correct saddle point in / will be recovered). The onset of a non-ergodic phase with many 
components (and thus the value of the dynamical transition temperature T c ), is signalled 
by the existence of a non-trivial limit for the order parameter, which is now the overlap 
of any two of these m copies in the carefully ordered limit lim m _>i lim 9 ^ + lirnzv->oo- (^ n 
a replica language what one needs to compute is the free energy within the 'one step rsb' 
Ansatz, where the n replicas are grouped into n/m groups of m replicas, an then expand 
the resulting free energy at n = around the point m = 1: F = Fq + (1 — m)F\. One then 
studies whether the piece F\ has a non trivial saddle point in the overlap.) 

For technical reasons which we shall not explain here, this criterion is known under the 
name of 'marginal stability'. It has been used in many instances for discontinuous spin- 



glasses. An insightful interpretation due to MonassonE3 shows how this procedure amounts 



to calculate a partition function restricted to a subset of configurations, chosen by a random 



pinning field. A related approach£23 uses another technique to calculate a partition function 
over a subset of configurations, where one starts from a typical configuration at given energy 
called the 'pivot', and calculates the free energy —1/(3 log(Z q ) associated to all configurations 
at overlap q from the 'pivot'. T c now appears as the temperature where this free-energy 
develops a local minimum for q 7^ 0. 

4 Glasses and spin-glasses without disorder 

As emphasized above, the Mode-Coupling equations which have been used with some suc- 
cess in the recent years to describe supercooled liquids 113 formally coincide with the exact 
equations describin g some mean-field spin-glasses, or the motion of a point particle in a 
random potential | 18 H 133 H 134 I (in large spatial dimension). This is a priori surprising in view 
of striking differences between the behaviour and the basic constitutive ingredients of spin- 
glasses on the one hand and structural glasses on the other hand. Let us mention a few 
rather obvious ones: 

- In spin-glasses, there exists some quenched disorder (e.g. the random position of the 
spin carrying atoms in spin-glasses, which do not move with time) which is absent in glass 
forming liquids, where the disorder is, in a sense clarified below, 'self-induced'. 

- The glass transition is a dynamical effect, basically defined as the temperature or den- 
sity at which the relaxation times reaches the experimentally accessible order of magnitude 
of hours. The existence of a true transition at lower temperatures is very controversial. In 
spin-glasses, on the other hand, there exists some (experimental L3t3 and numerical E3) evi- 
dence in favour of the existence of a second order phase transition (at least in zero magnetic 
field), with a power law divergence of the relaxation time and of the non-linear magnetic 
susceptibility. 

- Another important difference is the existence of a crystalline phase in structural 
glasses, which can in principle be reached in a very slow annealing procedure. This has no 

42 



counterpart in spin-glasses. 

In spite of these important differences, there has been recently an increasing convergence 
between the two fields, both at the theoretical level and experimentally, because of the 
existence of an aging regime at low temperatures. This chapter will review some of these 
points of convergence. 

4-1 Phenomenology of glasses: a few basic facts 

Many very different glass formers exhibit surprisingly similar properties when approaching 
the glass transition. The glass transition temperature itself is a purely conventional (and 
somewhat anthropomorphic) temperature where the relaxation time reaches a value of the 
order of 10 3 seconds. However a common experimental feature is the stretching and shoul- 
dering of the relaxation (of, e.g. the density fluctuations) as the temperature is decreased 
0'E3 towards T g . More precisely, the relaxation evolves from a simple Debye exponential at 
high temperatures (liquid) to a two-step process at lower temperature (supercooled liquid), 
where the correlation function first decays rather quickly to a 'plateau', and later departs 
from this plateau value on a much longer time scale t(T). Correspondingly, the frequency 
dependent susceptibility x"(u;) evolves from a one peak, high frequency, structure to a two- 
peak structure when the temperature is decreased. The second, low frequency peak (called 
the 'a-peak') shifts to lower and lower frequencies u; a = l/r(T) as T is lowered; the shape of 
the peak is furthermore strongly non-Debye, which reflects the fact that the time relaxation 
functions are non-exponential (and often fitted by stretched exponentials). The shape of 
the minimum lying in-between these two peaks has been the focus of an intense interest 
recently, essentially because one of the major predictions of mct is that, around a certain 
temperature T c , 

\"u) k j i — '" 



X"(Wmin) UrH W »^ 



^min j 



where a and b are two (positive) exponents related through Eq. (|54|). This behaviour 
reflects, in frequency space, the behaviour of the correlation function C(t) represented in 
Fig. (see Eq. (||)). 

The relaxation time r(T) grows extremely fast as the temperature is decreased, in 
general faster than the Arrhenius law exp(A/T). Systems for which t(T) are close to an 
Arrhenius behaviour are called 'strong', whereas systems for which the divergence is faster 
are called 'fragile' E3. For the latter systems, a widely used description of the experimental 
data (based on dielectric measurements, viscosity measurements, etc.) is the Vogel-Fulcher 
law 



r(T)~t e T - T o (79) 

which suggests that t(T) actually diverges when T — » To, i.e. that there is a true phase 
transition at T = Tq. This temperature furthermore appears to coincide with the tempera- 
ture (called Kauzman temperature) at which the extrapolated excess entropy (as compared 
to the crystal) would vanish. However, the divergence of t(T) is so rapid that t(T) becomes 
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larger than experimental time scales at the glass temperature T g which is often apprecia- 
bly larger than To. Is is thus difficult to claim the existence of a transition on the basis 
of the fit of t(T) only (see however l 136 l) In particular, othe r fu nctional forms, such as 
t(T) ~ £ exp((A/T) 2 ), also give reasonable fits of the data I1H1 - without invoking the 
existence of a critical temperature where t(T) would diverge 
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4-2 Discontinuous spin-glasses: a mean-field scenario for structural glasses 

In order to understand better why some theoretical ideas emerging from spin-glass mean- 
field theory might also be relevant for structural glasses, it is important to keep in mind the 
fact that, as we already emphasized in Section [|, there exist two different classes of mean- 
field spin- glasses, continuous and discontinuous. While the most conventional - continuous - 
ones provide a good starting point for the description of real spin-glasses, with a second-order 
phase transition where the Edwards- Anderson order parameter vanishes, some discontinuous 
models are characterized by a discontinuous static phase transition at a temperature T s , 
where the order parameter is finite just below T s . An extreme example of this type of 
spin-glass ordering is provided by the Random Energy model E3, which has zero entropy 
density in the whole low temperature phase. There are in fact many other such examples 
(see Section ([3, 2D), As discussed there, the discontinuous models generally possess a rather 
peculiar dynamical behaviour, with a dynamical transition temperature T c which is higher 
than the static one (and which coincides with the Mode-Coupling critical temperature). 
When approaching T c from above, the relaxation time t(T) diverges as a power law, but 
there is no singularity in the static thermodynamic quantities, which are analytic around 
T c (see Fig. ^). Thermodynamic singularities, including a jump in the specific heat, 
only occur at the lower temperature, T s . As explained above (Section 3.8), this behaviour 
originates from the fact that as soon as T < T c , the system starts aging for ever in a slow 
descent towards a state with a free-energy extensively higher than the one of the equilibrium 
state. Note that between T s and T c , the number of such metastable states is exponentially 
large, the associated configurational entropy (the 'complexity') being exactly equal to the 
free energy difference between the metastable states and the equilibrium (paramagnetic) 
state. 

Such a scenario can however only exist, strictly speaking, at the mean-field level, where 
nucleation barriers are infinite. In finite dimensional space, metastable states with a free 
energy density larger than that of the ground state have a finite lifetime, since the nucleation 
of large 'bubbles' of the ground state are always favoured. The gain in free energy for a 
bubble of size L scales as L , while the surface energy cost scales at most like L®^ 1 (if the 
interaction is short-ranged). Therefore the mean- field picture cannot survive as such in a 
finite dimensional system. However, one can argue that around the mean-field T c the bubble 
nucleation will be a slow process, leading to a rapid increase of the terminal relaxation 
time t(T), which only diverges when the static transition temperature T s is reached. This 
picture was proposed nearly a decade ago in an insightful series of papers by Kirkpatrick, 
Thirumalai and Wolynes Il3~0: within this hypothesis T c is a crossover temperature below 
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Figure 12: Relaxation time vs temperature in discontinuous spin-glasses. The right hand curve is the 
mean-field prediction, which gives a dynamical transition at a temperature T c above the static transition 
temperature T s . The left curve is a conjecture on the behaviour in finite dimensional systems: activated 
processes smear the dynamic transition. The relaxation time diverges only at the static temperature T a , but 
becomes experimentally large already around the glass temperature T g . 



which activated processes become important f] leading to a rapidly increasing (d la Voger- 
Fulcher 1123) relaxation time as the temperature is further reduced. The static ordering 
critical temperature T s appears as a Vogel-Fulcher or Kauzman temperature (below which 
freezing is complete), while the experimental glass temperature T g , where t(T) reaches 10 3 
seconds, lies somewhere between T s and T c (see Fig.12). 

This rather appealing idea however suffers from important theoretical loopholes. First of 
all, it relies on a model with quenched disorder, absent in structural glasses. As discussed in 
the next paragraph, this might not be too serious as this disorder might well be 'self-induced' 
by the system itself. One would actually like to be sure that the above nucleation arguments 
are indeed correct for finite dimensional versions of discontinuous models (there exist a few 
numerical simulations for the Potts glass, p-spin spin-glass and 'frustrated percolation' 
models IH!J~ti3) . The subtlety comes from the fact that the nucleation is in the present 
case rather peculiar, since the nucleating phase cannot be the ground state (otherwise the 
system would be completely frozen after a finite time, and would loose the contribution of 
the complexity to the free-energy), but rather another metastable phase with exactly the 



bubbles should grow 
] is not very clear to 



same free-energy density - which makes it hard to understand why t : 

at all. The meaning of the 'entropic driving force' invoked in Ref. : 

us, and, surprisingly, little progress has been made to support this conjecture t?J. Another 

picture, somehow related to that of Ref.^13], is that the complexity induces a microphase 

^Similar ideas can be found within the context of 'extended' MCT,liZI although the precise relation between 
the two pictures is not obvious to us. 
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Figure 13: The energy of the LABS model with open boundary conditions vs temperature. The top curves 
are the results from a Monte-Carlo simulation of a N — 401 spin system with logartihmically decreasing 
cooling rates (from Ref.[158]). The other curves are derived analytically with the fiduciary random system. 
T c is the dynamical transition temperature below which the fiduciary system freezes, and T s is the static 
transition temperature which cannot be found from a Monte Carlo, but is accessible in principle from exact 
enumerations of systems with a much smaller size. The lowest curve, corresponding to the 'crystalline' state, 
has been found only in the case of periodic boundary conditions and for some special values of N. 



separation into 'grains', each of a different metastable phase, with a certain temperature 
dependent size. The relaxation time would correspond to the time needed to a 'grain' to 
disappear or for a new grain to appear. In this sense, glass dynamics might have a lot in 



common with foams or microemulsions, as recently advocated in Ref.fM]. 



4-3 Self-Induced Quenched Disorder: Spin glasses without disorder 

An important obstacle if one wants to convert the above picture valid for some spin-glasses 
into a theory for structural glasses is the meaning of the qu enched disorder in the latter 
case. It turns out however that a series of recent works I113~E^3 has shown the existence 
of discontinuous spin-glass like behaviour in systems with frustration but without quenched 
disorder. These systems thus provide natural spin analogues of glass formers. Although their 
microscopic description does remain remote from that of structural glasses (in particular 
because they involve infinite range interactions), they provide at least an existence proof to 
the phenomenon of self-induced disorder, and their study is worth the effort, both for their 
intrinsic beauty and as a source of inspiration for modelling structural glasses. Furthermore, 
from an experimental point of view, Charge Density Wave systems (among others t5J) have 
recently been shown to behave very much like disordered systems B3, with however a very 
small density of impurities, suggesting that incommensurability effects alone (inducing some 
frustration) might be sufficient to generate 'self-induced disorder' 
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Monte Carlo simulations show that the low-temperature dynamics is highly non-trivial 
in all cases. Although an analytical study of the dynamical features could in principle be 
done along the lines described in Section 0, in many of the recent papers the shortcut was 
taken of calculating some dynamical properties using the 'pseudo-statical' approaches we 



described in Section 3.7. 



• Low autocorrelation binary sequences. 

We shall first concentrate on some spin systems with frustration but without disorder, 
which contain long-range interactions. These systems exhibit the same behaviour as the 
discontinuous mean-field spin-glasses, namely a dynamical transition at a temperature T c 
larger than the temperature T s of the static transition. The first example is the problem of 
'low autocorrelation binary sequences' (labs). This is an old and important pro blem from 
communication theoryll^B which was restated in physical terms by BernasconiE^H as follows: 
take a one dimensional chain of Ising spins <7j = ±1, i = 1...N. Compute the correlation 
function at distance k 

N 
C k = ^ a i a jSj,i+k (80) 

and define the energy function as 

-. JV-l 

Ei{a}) = W^T) J{ ° l ' (81) 

The interest in communication theory is to find the low energy configurations. In fact 
the ground state of this energy function does provide a sequence of bits which minimizes 
the two point correlations, and this is also useful for building a pseudo random number 
generator. Two versions of this problem have actually been studied, differing in the choice of 
boundary conditions. Due to the infinite range of the interactions, they present significative 
differences. A first version studied in Refs. ^a>££3>£i>3] has free boundary conditions, where 



the correlation function Ck is defined with the sum over the spin indices i going from 1 to 



N — k. Another choice, studied inlM3 is that of periodic boundary conditions, where one 
can define Ck through a sum over the spin indices i going from 1 to N — 1. As any other 
optimisation problem, this can be generalised to a finite temperature study by assigning to 
each sequence of spins a Boltzmann weight P({a}) = exp(—f3E({a})/Z. 

Let us first discuss the case with free boundary conditions. It was shown by Monte 
Carlo simulations that there exists a finite temperature freezing region in a temperature 
range around T ~ 0.1, with a weak cooling rate dependance of the low temperature energies 



(see Fig. 13). Computations spanning very long time scales have been perfomed at low 



temperature by using an efficient Monte Carlo algorithm l 158 l~H£9, and reveal a clear aging 



effect, characterized by a r/t w scaling (see Fig. 14). The smearing of the transition and the 
cooling rate dependence might be a finite N effect. In this case, the presence of 'traps' in 
phase space with a broad distribution of trapping timesll^3 is rather convincingly observed 
for finite N (see Fig.|l5|). 
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• A replica analysis for non-disordered problems. 

The analytical study of the LABS model is in itself very interesting. Despite its simplicity, 
we know of no direct solution. A rather indirect, but illuminating, way of proceeding is to 
replace the non-disordered LABS model by a 'fiduciary' model with quenched disorder. The 
basic idea consists in considering the model at hand as one special sample of an ensemble 
of systems containinipg quenched disorder. In the case with free boundary conditions this 



is achieved as follows l 145 H 147 1: one defines a 'disordered' correlation function 

C fc d = VM( fc V^-, (82) 



ij 



where M^ > is a matrix with random elements, equal to or 1, with the only constraint 

" . . M- (fc 



that ^2a M-a = N — k. The original problem is a particular choice of M^ k \ where the 



only nonzero elements are on the k diagonal. The hope is that this particular case is a 
generic case, and this is actually not at all obvious. (For example, it would be nonsense 
to claim that a ferromagnet is a special instance of a spin-glass with J^ = ±J couplings, 
where all Jij happen to be equal to +J: the ferromagnet is just a very atypical sample.) 
There is in fact quite a bit of educated guesswork involved in the choice of the ensemble 
of disordered system, of which the original model is argued to be a generic member - see 
below. In the present case, the original model is extremely frustrated due to the long-range 
and conflicting nature of the interactions, two features which are indeed retained by the 
Hamiltonian defined using Eq. (|82j). 

Now, the crucial remark is that if the model is indeed generic, its static properties can 
be obtained by means of the replica method, where the averaging is performed over the 
fictitious disorder. In the case at hand, the resulting free energy indeed turns out to be a 
good approximation of the original model in the high temperature, replica symmetric phase. 



This approximation actually corresponds to the one proposed by Golayll53 using different 
arguments; as seen from a high temperature expansion, this approximation is however not 
exact (but see next paragraph). Its main virtue is to predict the existence of a static phase 
transition at a temperature T s = 0.0476, below which a breaking of replica symmetry of the 
discontinuous type appears (see Fig. [b|). The low temperature phase is characterised by 
a residual entropy density which is linear in T, but small (less than 10 -5 per spin at T s ). 
From a glass point of view this phase transition can be seen as the resolution of an entropy 
crisis appearing at an extrapolated Kauzman temperature which is very close to T s . The 
prediction for the ground state energy density, Eq/N ~ 0.0202, is compatible with a large 
N extrapolation of the ground state energies found by exact enumeration on small samples 
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153 N < 48. On the other hand, it does substantially differ from the apparent ground 
state energy extracted from Monte Carlo simulations, even after extrapolating to very small 
cooling rate (see Fig.[l3|) . In fact, similar discrepancies have been seen and studied in detail 
before on several disordered spin-glasses (for instance, the binary perceptron problem 11^3) . 
This is again reminiscent of the above discussion of discontinuous spin-glasses and of the 
existence of a dynamical transition at T c > T s , where the Langevin (or Monte Carlo) 
dynamics gets trapped in metastable states with high energies -Ethres = E(T C ). In analogy 
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Figure 14: Two time correlation C(t w + r,t m ) obtained in an extensive Monte Carlo simulation of the LABS 
model of 400 spins with open boundary conditions (Ref. [158]), at a temperature T = .075, plotted versus 

r/tw The data exhibits a clear aging effect. 



with other discontinuous spin-glasses, one may thus expect that the dynamical transition 
where the energy freezes takes place at a temperature T c fixed by the marginality condition 



(see Section 3/7), i.e the temperature where a replica symmetry broken solution first appears, 
with a breakpoint (in Parisi's Ansatz) equal to m = 1. This leads to 1113 T c = 0.103, in 
reasonable agreement with the Monte-Carlo data (see Fig. |13| ). 



A model with a 'crystalline' state. 



In summary, we have shown that the LABS model with free boundary conditions provides 
an interesting example of a non disordered mean-field spin system, sharing many similarities 
with discontinuous spin glasses. However, the fiduciary Hamiltonian constructed using the 
disordered correlation Eq. (|8^) is only approximate, and does not, for example, give the 
exact free-energy in the high temperature phase. I t t urns out that the situation is under 
better control with periodic boundary conditionsl 147 H 148 l, In this case, the model furthermore 
exhibits, in some sense, a 'crystal' phase. 

The simplification comes from the fact that with periodic boundary conditions, the 
energy can be expressed as 

1 N 
E = J^E\^)\^ ( 83 ) 

q=l 



where we have introduced the Fourier series a(q) = J2j °j exp(2i7r<7j/./V). The problem 
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Figure 15: Energy per spin vs log 10 (time) in a single Monte Carlo run for a N — 400 LABS spin system at 

T = .075 (From Ref. [158]). 



remains non trivial since one has to remember that s(q) is constrained to be the Fourier 
transform of a binary sequence. In this case a clever choice of the fiduciary disordered 



system was found in till, which consists in substituting the usual Fourier components a(q) 



by a disordered version of the Fourier transform: 



(84) 



where U is a kind of random unitary matrix ]]. One can again compute the equilibrium 
properties of the disordered problem with the replica method. One finds, similarly to the 
case with free boundary conditions, t wo tran sition temperatures T s and T c . However, the 
periodic labs problem is in fact richer 
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- It can be shown that the disordered model has the same free energy as the original one 
to all orders in a high temperature expansion. The phase transitions, both dynamical and 
static, predicted by the replica solution of the disordered problem are in good agreement 
with the numerical simulations. The identity of the two problems in the low temperature 
phase, at the level of extensive thermodynamic quantitites, is however still a conjecture. 

- It turns out that the periodic LABS model possesses a very special ground state for 
some particular values of N. When N is prime, clever number theoretic properties can be 
used to generate a sequence of spins with a finite energy E, and therefore a vanishingly small 



There is a subtlety in this construction, namely the fact that one wishes to introduce a 'disordered 



Fourier transform' with a matrix U which is unitary, but also satisfies U- 



= U*j, where 



stands for 



complepc conjugation. The proper construction through random orthogonal transformation is described in 
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energy density E/N in the thermodynamic limit. (Actually for N prime of the form 4k + 3, 
k integer, the ground state energy is E = 1, and for N prime of the form TV = 4k + 1, 
is is found to be £ = 5). These special configurations are very difficult to find using 
Monte-Carlo like dynamics: the energy landscape has sometimes been described as a 'golf 
course' potential, in the sense that the ground state is an unexpected deep hole surrounded 
by rather unfavorable states. These ground states constitute the analogue of a crystalline 
state. If one simulates a LABS model with a 'magic N' (in the above sense), starting at zero 
temperature from the crystalline state, one finds that the low temperature specific heat is 
very small, and the energy has a discontinuous jump (first order transition) to the high 



temperature energy curve at a 'melting' temperature T m > T c (a sketch is given in Fig. 13). 
No such crystalline state has been found in the free boundary case c£J. 

• Speculations on the 'fiduciary' disordered Hamiltonian strategy. 

We have thus seen how some frustrated spin systems without disorder can be solved 
(approximately, or even exactly at least in the high T phase for the LABS with periodic 
boundary conditions), following a rather interesting strategy. This strategy consists in sub- 
stituting the original problem by a 'fiduciary' one with quenched disorder, and solving the 
disordered system using, e.g., the replica method to obtain the static properties and infor- 
mation about the transitions. There is unfortunately no systematic method of choosing the 
fiduciary model so far. The above two examples, or other models which have been stud- 



ied in a similar way lil3'E£iJ, show the importance of symmetry considerations in the choice 



of the fiduciary disordered problem, and suggest as a criterion that this disordered model 
should be as 'close' as possible to the original one in the high temperature (liquid) phase. 
This strategy is reminiscent of the very fruitful approach to energy levels in complex nuclei 
through the study of fiduciary random Hamiltonians with the proper symmetries Hill In 
our case we do not yet understand when such an approach may be successful or not, if it 
is only restricted to fin ite tirne d ynamics or if it does apply to thermodynamical proper- 
ties. In some cases (seel 133 H 148 l'E^ and below), the 'spin-glasses without disorder' explicitly 
involve pseudo random numbers in the sense that the spin couplings are deterministic, but 
very rapidly oscillating. This is much less obvious in the labs model, especially with free 
boundaries. Finally, as discussed in the next subsection, self-consistent (Mode-Coupling) 
approximations of non disordered models often lead to equations which are exact for some 
adequately chosen disordered systems. 

In a loose sense, one expects that the slow dynamics at low temperatures originates from 
some degrees of freedom which freeze and play the role of an effectively quenched disorder 
field for the other degree of freedom. The success of the present strategy might lie in the 
fact that the precise identification of these 'slow' variables is not necessary to understand 
the freezing transition in these models. 

• Other mean-field spin-glass models without disorder. 

There exist by now a few such examples of spin-glasses without disorder. Besides the 
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'fiduciary' Hamiltonian strategy, other approaches have been attempted, in particular some 
direct diagrammatic calculations - either static onescl3'E3 (basically through a resummation 
of the high temperature expansion) or using dynamical perturbation theory E2MLM Because 
of space limitation we cannot mention them all here. We should nevertheless point out the 
case of Josephson junction arrays b ecause it may allow a direct experimental investigation 



of many of the ideas discussed hereli^il'". This model involves two sets of N spins, a 



j 



exp(i(f)j) (living on the rows of a 2-D array of Josephson junctions) and a' k = exp(i<^.,) 
(living on the columns of the same array) which are normalised two component vectors 
(U(l) spins). These are coupled through the Hamiltonian 



j N 

H = -^Y. (exp&najk/N^oi + c.c.) (85) 

v 7V j,k=l 

where c.c. means complex conjugate. This can be realised by using a stack of two mutually 
perpendicular sets of N parallel wires. There is supposed to be a Josephson junction at each 
crossing of an horizontal wire with a vertical one. The variable (fij is the reference phase of 
the j th horizontal wire, while c/) k is the reference phase of the k th vertical wire. The system 
lies in an external transverse field H , which induces a phase shift per unit area a proportional 
to H. The high temperature expansion has been resummed for 1/N <C a <g 1, both for the 



static theorylilil and for the two-time correlation and response functionE£3, leading exactly 



to the same equations as those describing the p = 4 spherical spin-glass which we have 
discussed in Section ||. Again, the system undergoes two phase transitions, a dynamical 
one at a temperature T c larger than the static one T s . These dynamical equations are also 
equivalent to a particular case of the 'schematic' mode coupling equations of supercooled 
liquids. Josephson junction arrays may thus provide an interesting experimental playground 
to test directly the predictions of these theoretical studies. 

• Towards the description of the glass phase for interacting particles 

As is often the case in statistical physics, it is technically easier to study a phenomenon 
(here the glass transition in spin systems without disorder) on the magnetic case. A natural 
problem is to try to extend the kind of ideas that we have seen at work on these spin ana- 
logues towards more realistic problems where a glass phase appears without any quenched 
disorder, like glass forming liquids. 

One must thus describe interacting particles. The studies in this direction are still 
rather preliminary. Two routes have been recently explored. One is to consider interacting 



particles in large dimensional (D) spaces. In Ref. t£3] n point-like particles interacting 
repulsively are constrained to move on the vertices of a D-dimensional hypercube. The 
model can be mapped onto a modified 0(n) matrix model. It has a sharp glassy transition 
of the kind described previously. Amusingly, for n multiple of four, an unproven conjecture 
of Hadamard gives the 'crystalline' ground states. On the same line, similar properties are 



found in a modeled of hard spheres moving on the surface of a D-dimensional hypersphere. 

Another study retains the three dimensional nature of the problem, and starts from 

a reasonable analytic approximation of the liquid phase given by the Hypernetted Chain 

52 



Approximation (hnc). It turns out that one can generalize this HNC description through 
the introduction of a replicated theory as described in Section 3.7. This approachli^ yields 



a good analytical prediction for the glass transition temperature or density of hard or soft 
spheres systems. A precise description of the low temperature phase along the same lines 
is however still missing. 

4-4 p-spin models, Mode Coupling Theory of glasses, and its extension at low temperatures 

The above arguments suggest that discontinuous mean-field spin-glasses, in spite of the 
presence of quenched disorder, can provide a good starting mean-field theory for structural 
glasses. Among the most striking convergence between the two subjects, already alluded to 
many times above, is the equivalence (in the high temperature phase) between the Mode 
Coupling equations for glasses and the dynamical mean-field equations for p-spin-glass mod- 
els. This analogy was first noticed already long ago 113. It has been useful technically, in 
particular by transopsing the mathematical analysis of the Mode Coupling equations devel- 
oped for glasses IIZIEj to the stu dy of th e high temperature dynamics of p-spin systems c£! 



and manifolds in random media E23'tl3, e.g. the divergence of the relaxation time and the 
should ering q fth e relaxation. Following similar developments in the field of fully developed 
turbulence l 165 H 166 l, it has been realized more recently that the factorization property which 
is at the heart of the mode coupling approximation actuaLV becomes exact for certain sys- 



tems, which turn out either to contain quenched disorderly or some deterministic version 



of disordered systems, in the sense o f ha ving rapidly oscillating couplings which have sta- 



tistical properties of disordered ones £221. We shall first give a general flavour of why this 
is the case, and then turn to the implications of this general result to the low temperature 
extension of the Mode-Coupling equations, and its physical consequences. 

• Mode coupling approximation and hidden disorder. 

Both glass forming systems and turbulence can be described by some non-linear stochas- 
tic dynamical equation. In order to describe the essence of the mode coupling approximation 
and its relation with disordered systems, we shall explain it briefly on the simple case of a 
single scalar degree of freedom (j), with a Langevin dynamics 

^ = -M(*)*-f^ 3 + »7 (86) 

with initial condition <p(t = 0) = 0. The thermal noise rj is a Gaussian noise rj with 
{vify ) = an d (??(*) ??(*') ) = 2T 5(t — t') (in the following the brackets will always denote 
an average over the realisations of the Gaussian white noise rj). The coupling constant g 
serves as a book-keeping parameter to set up a perturbative expansion. This expansion can 
either be well-behaved or ill-behaved depending - say - on the dimension of space. It is in 
any case rather useless when g = 0(1) if it cannot be resummed in one way or another. The 
mode coupling consists in a 'one-loop' self-consistent perturbation theory. This amounts to 
resumming a particular (infinite) set of terms in the perturbation expansion. In this way, 
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non-trivial self-consistent equations are obtained, which enable one to peep into the strong 
coupling regime, through an approximation which is however not easily controlled. 

Setting Rq = [(J,(t) + ^] _1 , which gives Ro(t,t') = exp ( — J t , dr h(t)), the perturbative 
expansion for </>(£) is easily written as: 

<p(t) = R <8 r] - jL R ® {R & rj • R ® rj . R <g, v } + ... (87) 

where (8) means a time convolution: (Rq (8 /)(£) = / dt'Ro(t,t')f(t') and • is a simple 
product. 

The lowest non trivial order perturbative expansion for the correlation function C(£, £') 
and the response function R(t, £') is easily written in terms of the kernels X(£, £') and D(t, £') 
through the Dyson equations: 

R(t,t') = Ro(t,t f ) + / dti I * dt 2 Ro(t,h) £(£i,£ 2 ) R(t2,t') , (88) 

C(£,0 = [ dh [ dt 2 R(t,t 1 )D(t 1 ,t 2 )R(t',t 2 ). (89) 

jo jo 



The mode coupling approximation for this problem amounts to an approximation of the 
kernels £(£,£') and D(t,t') where one takes their values at order g 2 and substitutes in them 
the bare propagator Rq and the bare correlation by their renormalised values. This gives 
the following self-consistent equations: 

£(£,£') = ^-C 2 (t,t')R(t,t>) 

D(t,t') = 2T6(t-t') + ^[C(t,t')} 3 , (90) 

This approximation neglects 'vertex renormalisation'. 

The problem is of course to try to control this procedure. An important step in this 
direction is to identify a model for which the self-consistent equations are exact. The ba- 
sic remark (first made by Kraichnan in the context of turbulenceli^J, where the analogous 
method is named direct interaction approximation) is that the diagrams retained by the 
mode coupling approximation are precisely those which survive if one considers the fol- 
lowing disordered problem. First, one upgrades <f) to an N— 'colour' object (f) a , where 
a = {1,2,..., JV}. The equation of motion Eq. fl86D is then generalized to: 



— — = -fx(t) 4> a -ig 2^ J al 3~/5 4>b4>"/4>6 + Va (91) 

m (3<J<8 

with independent noises r\ a . The couplings J a /3yS are quenched, i.e. time- independent, 
Gaussian random variables of zero mean and variance J^s-yS = V-^ 3 - in the large N 
limit, the correlation: C(t,t') = -^ J2 a =i (4>a(t)<fta{t')} (where the overline denotes the 
average over the random couplings J a /3-ys) and the response: R(t , £') = jj J) a=1 ( Q^ a L\ )\. n =o 



precisely obey the mode coupling approximation equations I12j, Eqs. (|8g| ) and (| 
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A particle in a random potential. 



The same construction can be generalised 11211 to an arbitrary nonlinearity F{4>) substi 



tuting the Jy <ft s in Eq. fl86|) . Therefore one finds that the general schematic mode coupling 
equations developed in the study of glass forming liquids can be derived exactly from the 
Langevin dynamics of N continuous spins (f) a , of the type: 

-M(*)0a-^Sr^+»7«. (92) 



Of 



This disordered multispin Hamiltonian is precisely the generic mean-field problem ( j44| ) 
which we studied in Sect. ||, which can also be seen as describing a particle evolving in an 
N — ► oo dimensional space in a quenched random potential V[{0}]. 

This interpretation is rather appealing. Let us introduce the following highly simplified 
picture of a glass: the motion of a given particle can be thought of as taking place in a 
random potential created by its neighbours. Since the motion of the molecules is extremely 
slow at low temperatures, one can assume that this random potential has a static component, 
in the spirit of the 'self-induced quenched disorder' scenario which we discussed above 11^3. 
In large dimension of space, one can establish the exact equations relating the two-time 
correlation function C(t w +r, t w ) = (f(t w +T)-f(t w )} (where r(t) is the position of the particle 
at time t), and the two-time response to an external force R(t w + r,t w ). For temperatures 
higher than T c , we have seen that C and R are actually tti, and furthermore that the fdt 
R(t) = — ^0(t) q-J' is obeyed. As noted in Section 3.2, one can then eliminate R(t) and 
find an equation for C{t) which is precisely the schematic Mode-Coupling equation, with a 
kernel related to the correlation of the random potential. 

Hence, the physical content of the (schematic) mct is clear: it is a mean-field description 
of a single point in a static quenched random potential. The important point is thus that 
mct implicitly assumes the presence of some quenched disorder which should rather, as 
discussed above, be 'self-induced' by the dynamics itself. In a sense it looks rather similar 
to the introduction of fiduciary models discussed before. 

• Mode coupling at low temperatures 

Coming back to the general equations relating C and R, one can now postulate that 
they are the correct generalisation of the schematic Mode Coupling equations for two time 
quantities, and investigate the 'glass' phase T < T c . The results of Section 3 are thus directly 
applicable. In particular, the correlation and response function cease to be functions of r 
only. More precisely, C(t w + T,t w ) can be written as the sum of a stationary contribution 
C ST (r) which only depends on r, and an aging part which depends on the ratio t w /(t w +r), 



(or a generalisation thereof, see Section 3.6): C AG (t w /(t w + r)). The expected shape of the 
correlation function in the glass phase is thus given in Fig. |I| 

The same decomposition holds for the response function, and the aging parts of C and 
R are related by an 'anomalous' fdt, where T is replaced by an effective temperature T/X, 
with X < 1. 
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In more physical terms, this means that for a finite waiting time t w after the quench 
below T c , one expects that the susceptibility xi^^w) still exhibits two peaks: a high fre- 
quency /?-peak very similar to the high temperature (T > T c ) one, and a low frequency 
a-peak which reaches a maximum at a frequency w Q ~ l/t«j, which thus progressively dis- 
appears as t w — ► oo. An interesting prediction of this low temperature extension of mct 
is that the high frequency part of the aging a-peak behaves as (uit w ) , while the low fre- 
quency 'foot' of the /3-peak behaves as uj a , with the following relation between a, b, and 
X (see (H)): AT 2 [1 + b]/T[l + 26] = T 2 [l - a]/T[l - 2a]. This equation generalizes the 
well known mct relation (|54]) between a and b for T > T c , for which X = 1. It would be 
extremely interesting to try to test these predictions experimentally, taking care of the fact 
that the above picture is only valid insofar as t w is small compared to the relaxation time 
t(T), such that 'activated' effects might indeed be neglected. 

5 Conclusion. Where do we stand ? 

Let us now summarize some of the most important ideas developed in this review and 
discuss some open problems. 

We have tried to show that aging effects are not spurious, irreproducible artifacts of 
non equilibrium situations, but rather an unescapable feature of systems characterized by a 
very large relaxation time, because time-translational invariance breaks down and the well 
known fluctuation-dissipation theorem has to be modified in a non trivial way. However, 
the asymptotic aging regime where all times are large reveals some universal properties; in 
particular, the 'effective' relaxation time becomes a time dependent notion and grows with 
the waiting time. The detailed investigation of these aging effects actually offer a unique 
way to probe the phase-space structure of complex systems. Prom that point of view, a 
system in equilibrium is 'dead'. 

We have described several simple models where aging can be described in detail: coars- 
ening models, 'trap' models and mean-field spin-glass models: (which turn out to give 
dynamical equations in exact correspondance with the 'Mode-Coupling' description of su- 
percooled liquids). Coarsening leads to aging in the correlation function but not in the 
response function: the fdt is most strongly violated in the aging regime. This scenario 
can thus only include aging in the response of spin-glasses or polymer glasses as a transient 
effect even in an infinite system - quite apart from the fact that it is not obvious what 
would actually grow with time in these systems. 

Mean-field spin- glasses provide an interesting testing field. One generically finds aging 
in the low temperature phase of these models, although two very different categories of 
systems emerge. In the discontinuous case (corresponding in the high temperature phase 
to 'model B' of mct), a dynamical temperature transition is found above the equilibrium 
transition. Throughout the low-temperature phase the asymptotic energy-density arrived 
at after any type of cooling is higher than the equilibrium energy-density The two-time 
plane breaks up into two sectors, which correspond to the stationary and aging dynamical 
regimes. On the contrary, for the continuous case (such as the Sherrington-Kirkpatrick 
model) static and dynamical transition temperatures coincide, as do the asymptotic out 
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of equilibrium energy density and the equilibrium one. The two-time behaviour is much 
more complicated, reflecting in some way the subtleties of Parisi's ultrametric equilibrium 
solution. 

Interestingly, if one views mean-field glassy dynamics as the dynamics of a point in 
a many dimensional rugged phase-space, one finds that the basic mechanism for aging is 
germain to that of simple coarsening: because of high dimensionality of phase-space the 
system starts near the border of a basin of attraction, and remains forever undecided about 
where to go. 

It is of course crucial to know how these mean-field predictions are modified when one 
goes beyond mean-field and studies finite dimensional systems. From a theoretical point 
of view, the difficulty comes from the fact that activated processes, which are effectively 
absent in mean-field, come into play when the dimension is finite. For example, the in- 
finitely long-lived metastable states in mean-field acquire a finite relaxation time in finite 
dimension, through bubble nucleation. The dynamical transition, which corresponds in the 
language of supercooled liquids to the Mode-Coupling critical temperature, is thus smeared 
out. Similarly, a particle in a random potential in finite dimension reaches a local minimum 
after a finite time, beyond which thermal activation starts playing an important role; con- 
tributions to aging of a somewhat different nature — such as those described by the 'trap' 
picture -- then set in. The relation of these trap models to the aging dynamics of real 
spin-glasses or other disordered systems such as pinned vortex lines, dislocations, domain 
walls, or polymers is however still rather tentative. 

For the same reason (absence of activated effects), mean-field models are not suited to 
describe cooling rate dependent effects, which can be very large in some disordered systems 
and in glass forming liquids. 

Therefore, new theoretical ideas which would allow one to extend the previous ap- 
proaches to finite dimensions (bearing in mind that activated effects cannot, in general, be 
accounted for within perturbative schemes) are clearly desirable. Returning to the question 
of aging in the response, the FDT-violation factor X in finite dimension is particularly in- 
teresting: this factor — that is related to effective temperatures in the system — is found 
to be non-trivial (0 < X < 1) in the disordered mean-field models; what is the situation in 
finite dimensions? Some numerical resultsc3'E3 suggest that X does indeed remain non triv- 
ial in finite dimensions, while arguments based on coarsening pictures would suggest that 
asymptotically, X = 0. This might actually be a clear-cut dynamical distinction between 
droplet like coarsening pictures and mean-field like pictures of real spin-glasses, which surely 
deserves more efforts - as stated above, spin-glasses, after all, do exhibit aging in response 
functions. 

We have also discussed rather at length the relation between glasses and disordered 
systems, noting that: 

- some models without disorder behave very much like disordered systems and can 
actually be theoretically described as such, and 

- some approximation schemes for systems without disorder, but strongly interacting, 
lead to dynamical equations (i.e. the Mode-Coupling equations) which correspond to the 
hidden assumption of the existence of some quenched disorder and are actually exactly the 
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equations describing disordered mean-field spin-glasses. 

Correspondingly, many of the ideas developed to describe aging in disordered systems 
are de facto also relevant for glasses. We have mentioned in particular how aging should 
manifest itself, within the Mode-Coupling scenario and for times smaller than the relaxation 
time as a waiting time dependent a-peak. A detailed analysis of these aging regimes should 
enable one to distinguish, again, between mean-field like descriptions and activated, 'trap' 
like, models; or perhaps understand how both mechanisms are blended. 

Several questions still remain completely open. In particular, is there a general criterion 
allowing one to understand when a 'complicated' system can be described as disordered ? Is 
this description only viable at finite times, where the specificity of the system at hand has 
not yet had time to reveal itself ? A possibility is that this time scale can only be infinite for 
mean-field like models, such as the labs model or the long-range Josephson array. Related 
to this is of course the lurking question of the very existence of a true glass transition in 
any short range system without disorder. 

The question of the existence of a true phase transition is however, after all, not crucial 
to understand the out-of-equilibrium properties, as it is displayed in so many physical 
systems. In this respect, it is plausible that mean-field models provide in general a much 
better starting point for the finite time, dynamical properties of real systems, than for its 
long time, equilibrium properties. 
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